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f^ ' Abstract. We study the surface tension and the phenomenon of phase coexistence for the Ising 

/— ^ ^ model on Z'' (d ^ 2) with ferromagnetic but random couphngs. We prove the convergence in 

^», 1 , probability (with respect to random couplings) of surface tension and analyze its large deviations 

: upper deviations occur at volume order while lower deviations occur at surface order. We study 

the asymptotics of surface tension at low temperatures and relate the quenched value r' of surface 

^ ^ tension to maximal flows (flrst passage times if d = 2). For a broad class of distributions of the 

■^U( , couplings we show that the inequality r" ^ r' — where t'^ is the surface tension under the 

averaged Gibbs measure - is strict at low temperatures. We also describe the phenomenon of 

phase coexistence in the dilute Ising model and discuss some of the consequences of the media 

randomness. All of our results hold as well for the dilute Potts and random cluster models. 
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A considerable amount of work permitted to understand on a rigorous basis the phenomenon 
of phase coexistence in models of statistical mechanics like the Ising model. Phase coexistence 
in the Ising model was first described in the pioneer work [20] . in the two dimensional case and 
at low temperatures. The construction was then simplified |37j and extended up to the critical 
temperature [551 [551 [SO] j still in the two dimensional case. The generalization to higher dimensions 
was achieved later thanks to the L^-approach [5j^ Jjj. The interested reader will find pedagogical 
presentations of the problem and the methods in the course [TU] and the review [S] . 

The present work is concerned with the phenomenon of phase coexistence for the dilute model 
with random (ferromagnetic) couplings. The random couplings model either rare defects in the 
media, either intrinsic randomness . As an example, quenched alloys made of magnetic materials 
have intrinsic randomness since the strength of the interaction between two spins depends on the 
nature of the two corresponding atoms. 

In order to describe rigorously the phenomenon of phase coexistence in presence of phase 
coexistence, we followed the same plan as in the above mentioned works. In a first step we 
established a coarse graining for the model [44] . In a second step - the present one - we study 
surface tension. The combination of these tools allow us describe the phenomenon of phase 
coexistence in the presence of random media. 

Before we turn to the presentation of the model and of our results, we would like to stress two 
consequences of the media randomness on the phenomenon of phase coexistence: first, it is the 
case that the shape of crystals are smoother than in presence of uniform couplings. Second, we 
give an insight to the expected localization phenomenon of the crystal which is determined by the 
realization of the media under averaged Gibbs measure. 

The organization of the paper is as follows. In Section [T] below we introduce the model and give 
a complete summary of our results on surface tension, its low temperature asymptotics (maximal 
flows) and phase coexistence. Proofs and intermediate results are given in the three corresponding 
Sections [1 [3] and H 

1. The model and our main results 

1.1. The dilute Ising model. The canonical vectors of M.'^ are denoted (6^)^=1...^ and for any 
X = X]r=i ^i^i ~ (^1' ■ ■ • ' ^d) & I^*^ we consider the following norms on R"^: 

d / d \ 1/2 

(1.1) ||a;||i = ^ |xi|, ||a;||2= ^a^n and ||a;||oo = max |a;j|. 

Given x,y € Z"^ we say that x,y are nearest neighbors (which we denote a; ~ y) if they are at 
Euclidean distance 1, i.e. if ||a; — y||2 = 1- To any domain A C Z'^ we associate the edge sets 

(1.2) £:(A) = {{x,y}:x,yeA and x^y} 

(1.3) and ^""(A) == {{x,y} : x e A,y e Z'^ and a; - y} . 

We consider in this paper the dilute Ising model on Z'' for cf ^ 2. It is defined in two steps : 
first, the couplings between adjacent spins are represented by a random sequence J = {Je)eeE{Z'^) 
of law P, such that the (^e)e6_E(z<') are independent, identically distributed in [0, 1] under P. Then, 
given A C Z*^ a finite domain and a spin configuration cr G Sj, where 

Sj = {ct : Z'' ^ {±1} : a, = l,Vz ^ A} , 
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we let 

(1.4) H-l^+{o)^~ J2 ^-''-''y 

e={x,y}eE^'{A) 

the Hamiltonian with plus boundary condition on A. The dilute Ising model on A with plus 
boundary condition, given a realization J of the couplings, is the probability measure ^^ on Sj 
that satisfies 

(1.5) 4+({^})= 1 expf-|i/;(-+(a)), VaeS+ 

^a',/3 \ ^ / 

where (3 ^ is the inverse temperature and Z'^ „ is the partition function 
(1-6) z'^+= V e^pf-^Hi+ia] 



Consider 






J,+ 



(1.7) mp = lim E/i. (ctq) 

the magnetization in the thermodynamic limit, where A^v is the symmetric box A^v — {—N, . . . , N}''' 
and E the expectation associated with P. When mp > the boundary condition has an influence 
on the spins at an arbitrary distance. In the region to^ > we say that the Ising model has 
two phases because the structure of the spins under E^^'"*" (the plus phase) is not the same as 
the structure of the spins under E^^'^ (the minus phase), where /i^'~ corresponds to the minus 
boundary condition. 

It is shown in 1 that the dilute Ising model undergoes a phase transition at low temperature 
when the random interactions percolate. In our settings, this means that the critical inverse 
temperature 

(1.8) Pc = ini{P^O:m,}>0}, 

which is never smaller than (3^™° - the critical inverse temperature for the pure Ising model ( J = 1) 

- is finite if and only if P( Je > 0) > Pc{d) where Pc{d) is the threshold for bond percolation on Z"^. 

The aim of the paper is to understand the mechanism of phase coexistence in the dilute Ising 
model, hence we will consider in the following a distribution F of the couplings such that P( Je > 
0) > Pc{d) and an inverse temperature /3 > /3c- However, some of our results hold on a possibly 
stronger assumption P > (3c ^ (3c where (3c is the critical inverse temperature for slab percolation 

- see (|1.14p below - as this assumption allows us to use the renormalization framework of f3? . 

1.2. The Fortuin-Kasteleyn representation. The study of surface tension for the dilute Ising 
model will be led under the random-cluster model that corresponds to the measure /i^' ^ • We call 

n^ {oj : E{Z'') ^ {0,1}} 

the set of cluster configurations on E{Z'^), and for any oj G il and E C E{Z'^) we call lu^e the 
restriction of w to E, defined by 

('^l^^^ - \ else. 

The set of cluster configurations on E is VLe = {lo\Et^ G ^}- Given a parameter q ^ 1 and an 
inverse temperature /3 ^ 0, a realization of the random couplings J : E{Z'^) -^ [0, 1], a finite edge 
set E C E{'L'^) and a boundary condition tt G fi^c we consider the random cluster model ^e^o^ 
on ^E defined by 

(1.9) cj,^-^^({^})= 1 J|p-=(l_p^)i-=xg^5^(-), V^Gf7^ 
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where Pe = 1 — exp(— /3Je), C^(w) is the number of clusters of the set of vertices in 1/ attained 
by E under the wiring w V tt such that (lo V 7r)e — niax(ti;e, TTg), and Z^^n^ is the renormahzation 
constant making ^ ^^ a probabihty measure. 

For convenience we use the same notation for the probabihty measure ^^j^ and for its 
expectation. Most often we will take either tt = /, where / is the free boundary condition : 
/e = 0,Ve G i?"^, or TT = ui where w is the wired boundary condition : Wg = l,Ve G E"^ . When 
the parameters q and (5 are clear from the context we omit them. Given 7?, a compact subset 
of W^ (usually a rectangular parallelepiped) we denote by ^.^ the measure ^p!U „d. on the 

cluster configurations on £'(7?. n Z''), where 7?, stands for the interior of TZ. In particular, for any 
(7, /i : 17 -^ R the quantities ^-^ [g) and ^^ {K) are independent under P. 

The connection between the dilute Ising model A'^'a and the random-cluster model was made 
explicit in '22]. Consider the joint probability measure 



vl/^;+({(a,^)}) = ^^ n (Per(l-Pe)'-"% V(a,a.)GE+xf7, 



where pe = 1 — exp(— /3Je), a < lj is the event that a and w are compatible, namely that oj, 
7y,Ve = {x,2/}G£;"'(A),andZ;(;+ 

i. The marginal of '^ ^ a on the variable cr is the Ising model /i^' 

^B(A),/3 



1 =^ cTrr = crj,,Ve = {a;,j/} G £'"'(A), and Zj^g is the corresponding normalizing factor. Then, 

ii. Its marginal on the variable w is the random-cluster model 'I'^^'-, /^ with wired boundary 

condition w and parameter q — 2. 
iii. Conditionally on to, the spin a of each connected component of A for uj (now cluster) is 

constant, and equal to -|-1 if the cluster is connected to A'^. The spin of all clusters not 

touching A"^ are independent and equal to +1 with a probability 1/2. 
iv. Conditionally on ct, the edges are open (i.e. LUe ^ I for e ~ {x, y}) independently, with 

respective probabilities PeSa-^,a-y- 
According to point ii and iii we can study surface tension for the Ising model under the Fortuin- 
Kasteleyn representation. This representation allows to study at the same time the surface tension 
and the phenomenon of phase coexistence for the dilute Ising model {q = 2), but also for dilute 
percolation {q — 1) and for the dilute Potts model {q G {3,4, . . .}). 

An important benefit of the representation is that it makes possible the use of the comparison 
inequalities for the random cluster model. We say that a function / : He -^ K^ is increasing if, 
for all LUjiv' G O^; one has to ^o i-^' ^ f{^) ^ /(^') where ^ji stands for the product order on ^Ie- 
It was shown in |2^ that: 

i. For any h : Qe —^ R^ increasing, ^E^n'^i^) ^^ ^ non-decreasing function of J, /? and tt. 
ii. (FKG inequality) For any g, h : He -^ K^ increasing, 

(1-10) ^iyiah) ^ ^■'Eyig)<^'Ey{h). 

iii. (DLR Equation) For any E' C E and uj' G VLe' , 

(1-11) *^r/(-h^'=-')=*MTi''- 

Finally, let us recall the assumption of slab percolation, that is the basis for a renormahzation 
framework in the dilute Ising model (44j- When d ^ 3, we say that slab percolation occurs under 
£$'■'■« if, for large enough H, 

(1.12) inf inf E^-'/ (x <^ y") > 

where Sl,h is the slab Sl,h — {1, • • ■ , L}'^"^ x {1, . . . H}. When d = 2, we say that slab percolation 
occurs when there exists k : N* i-^ N* with limTv^oo k{N)/N = such that 

(1.13) lim E$'q (there is an horizontal crossing for uj) > 0. 
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The critical inverse temperature for slab percolation is 

(1.14) (3c = inf <j3^Q: slab percolation occurs under E<i>^'''''^ > , 

it satisfies (3c ^ (3c where (3c is the critical inverse temperature for phase transition in the dilute 
Ising (resp. Potts) model. We believe that (3c and (3c do coincide. Upper bounds on /3c are derived 
in [3J from the argument of [1]. The technical assumption (3 > $c allows us to use a coarse 
graining, which is a fundamental tool at the moment of defining the local phase of the dilute Ising 
model (see Theorem 5.7 in [Sj, or Section [3] below) . 

1.3. Surface tension. One of the main issue we address in this paper is the behavior of surface 
tension and the influence of the random couplings. We consider the surface tension in large 
rectangular parallelepiped oriented along some direction n e S'^~^, where S'^~^ is the set of unit 
vectors of M.'^. The other axes of the parallelepiped are represented by 5 G S„, where 

§„ = < \^[±l/2]Mfe; (mi, . . . , Mfe_i, n) is an orthonormal basis of M'' > 

U=i J 

is the set of d — 1 dimensional hypercubes of side- length 1, centered at 0, orthogonal to n G S"^^^. 
Finally, we call a; e M'' the center of the rectangular parallelepiped and L, H its side-lengths, and 
denote finally by 

(1.15) Tl^^L,H{S.ri) = x + LS+[-H,H]n 

the rectangular parallelepiped centered at x, with basis x A- LS and extension 2H in the direction 
n (see Figure [T]). The discrete version oi TZis TZ — TZ C] Z"^ and the inner discrete boundary of TZ is 



dn= lyen:3zeZ'^\n,zr^y\ 



For any TZ as in (|1.15|) we decompose dTZ into its upper and lower parts d^TZ = {y G dTZ 
iy-x)-n^O} and 9"7^ ^ {y e dTZ : {y - x) ■ n < 0}. 



Figure 1. The rectangular parallelepiped 7^a;^L,i/ (5,71). 

In the context of statistical physics, the surface tension is the excess free energy per surface unit 
due to the presence of an interface. The surface tension in TZ thus quantifies the probability of 
observing the plus phase in the upper part of 7?, and the minus phase in the opposite part under 
the measure /i^ with free boundary condition. It is more convenient to formulate the definition 
under the random cluster model, where we translate the event of phase coexistence into an event 
of disconnection. 

Definition 1.1. Let 7?, be a rectangular parallelepiped as in p.l5p . The event of disconnection 
between the upper and lower parts of &TZ is 

(1.16) Pk = (w G fi : 9+7^ <% 5"7t| 
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and the surface tension in TZ is 

(1-17) r^ = --^^log$^-(I?^), 



^J — i^„.,ft-^'"' 



We denote by J™'" and J'^'^^ the lowest and largest values of the couplings according to the 
support of P, that is to say : 

J""" = inf{A ^ ; P( Je < A) > 0} 
and J™^'' = sup{A ^ : P(Je > A) > 0}. 

We also denote by r^"" (resp. t^'^'^) the value of the surface tension in TZ corresponding to the 
constant couplings J = J""" (resp. J = J'^^^). We have: 

Proposition 1.2. Let TZ be a rectangular parallelepiped as in il.l5\} . with L,H ^ 2yd. The 
surface tension r^ is a non- decreasing function of J and (3. It is a non- increasing function of H . 
With probability one under P, 

where Cd < oo depends on d only. 

As in the uniform case pSl, surface tension is sub-additive (see Theorem 12. II) . and this implies 
convergence in probability of r.^(n). 

Theorem 1.3. There exists Tp{n) ^ 0, the quenched surface tension, such that, for all P ^ and 
n e S'^-\ 

N^^oo ^^o,«,5«(5,n) = ^^(") «^ ^-probability 
whatever is 5 G S^ and S > 0. 

Similarly, the surface tension for the constant couplings J™™ and J™^'^ also converge and we 
denote by T'^"'^{n) and T'°^'^^{n) their respective limits. 

The sub-additivity is of much help for controlling the order of deviations from the quenched 
value of surface tension rl (n) . Upper large deviations happen at a volume order hence they have 
no influence on the phenomena we study here: 

Theorem 1.4. For any e > and 6 > 0, 

We will be more concerned with lower large deviations. These are possible when T™^'^{n) < 
T'^{n). The inequality is known to be strict only in two specific cases: when J™'" = and /3 > /3c, 
it is the case that T'^"^{n) < T'^{n) because T™^'^{n) = 0, while the coarse graining [33] implies: 

Proposition 1.5. Assume (3 > Pc- For any n e S'^~^ , T'^{n) > 0. 

When P(Je > J™'") > Pc{d) and f3 is large enough, the strict inequality T'^"^{n) < r'^iri) also 
holds, cf. CoroUarv 1 1 . 141 below. 

When lower large deviations occur, they have at most surface order. It is another consequence 
of sub-additivity that: 

Theorem 1.6. For every n e S'^^^ and f3 ^ 0, t > r'"'"(n), the limit 
(1-18) /„(r) =lim-^logP(r4„^„(5,„) ^r) 

exists in [0, -t-oo) and does not depend on S > 0, nor on S £ §ri- In is continuous, convex non- 
increasing, and In (r) = for t ^ t'^ (n) . 
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Figure 2 . Large deviations of surface tension. 



(1.19) 



For convenience, we extend the definition of /„ letting 

f oo for T < r™'" 



Ir.{r) 



lim£_o+ In {t + e) at r = r" 



(n) 
(n). 



In order to show that lower deviations are exactly of surface order, we need to prove that In 
is positive on the left of r'. We developed an argument based on measure concentration coupled 
with a control of the length of the interface. In the case of the Ising model at low temperatures 
we could establish a first control: 

Theorem 1.7. Assume q ~ 2 and J™'" > 0. Then, for j3 large enough there exists c > such 
that, for all r > Q, 

(1.20) Ir,{T''{n)-r)^cr''. 

In the general case a careful adaptation of the method yields: 

Theorem 1.8. For every n £ S'^^^ , for Lehesgue- almost all /? ^ 0, 



(1.21) 



If3AT-l{n)-r) 
limsup J > 0. 



These quadratic lower bounds on the rate function generalize common controls for directed 
polymers models jl6[ |B], which were introduced in order to represent interfaces in the two- 
dimensional Ising model with random couplings at low temperatures [271 I16j . Its is probable 
however that the quadratic order in the former Theorems is not optimal in two dimensions, as the 
comparison with directed polymers suggests that 



(1.22) 



In{T'^{n)-r) ^ cr^/^ and ^ = -. 



This scaling has been established rigorously for the zero-temperature limit of directed polymers: 
last passage percolation, for a geometric distribution of the passage times, see Theorem 1.1 and 
(2.23) in [31]. 

Some of our results on phase coexistence and on the dynamics of the dilute Ising model [45j 
require that lower large deviations are actually of surface order. An easy but important consequence 
of Theorem ll.Sl is: 

Corollary 1.9. The lower large deviations are of surface order when j3 i— > T?(n) is left continuous. 
Hence the set 

(1.23) A/} = |/3 ^ : 3n e 5^"^ 

is at most countable. 



and 7' > such that Ip^ri{Ti{n) — ?-) = > 
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We end the presentation on surface tension with the definition of the surface tension under the 
averaged Gibbs measure. It is the Fenchel-Legendre transform of /„, 

(1.24) T^(n)= inf {At + /„(t)} 

which coincides with the surface tension under an average of Gibbs measures: 



1 



(1.25) .^(n) = ^linr^-__logE( $-» (P, 



.J.w 



^TZ" 



A>0,neS"^^i 



where TZ^ = TZo^N,5N{S,n), as shown in Proposition 12.21 The particular case A = 1 corresponds 
to the usual notion of surface tension under the averaged measure (or annealed surface tension) 
and we denote T°'{n) = T'*'=^(n). 

The asymptotics of r'^/A as A ^ or +oo are given in Proposition [531 An important question 
about the surface tension under the averaged Gibbs measure is whether the random media is able 
to turn Jensen's inequality 

(1.26) r^(n)^AT9(n) 

into a strict inequality. A partial answer to this question is given in the next Section. Let us 
explicit the connection between the strict inequality T''^{n) ^ AT'^(n) and the asymptotics of /„ on 
the left of T'^{n): the opposite of the slope of /„ on the left of r'(n) is exactly 

(1.27) a„ = sup{A> 0:T^(n) = Ar«(n)}, 

with the convention that sup0 — 0. 

Finally, as in the non-random case, the homogeneous extension of each of these notions of surface 
tension r'^, r^, t™™ and r™'*'' are convex and continuous (Proposition [ 



1.4. Low temperature asymptotics. The low temperature asymptotics of surface tension 
permit to give a more precise insight into the properties of surface tension in random media. First 
we need to introduce the concept of maximal flow through the capacities J, where J = (Je)ee_E(z<i) 
is the family of random couplings introduced in the former section. Here we give only a brief 
overview of maximal flows. The reader is invited to consult [331 [31] for a pedagogical introduction. 
Recent results on maximal flows, including large deviations, can be found in [i ^ [3 ^ IT7]. 

We will use an analogy for describing maximal flows. Imagine a liquid which has to cross a 
lattice made of tubes with limited capacity. Then, the maximal flow, in a given direction, is the 
quantity of liquid that can flow through the lattice, per unit of surface. 

Given a rectangular parallelepiped TZ as in p.lSp and I C EilZ), we consider the event that w 
is closed on /: 

Z7 = K = 0,Vee/}. 
We say that / is an interface for TZ \i Zi d Du and Ve e /,Z/\{e} /^T^tz- In other words, / is 
an interface for TZ if the disconnection on / is enough for disconnecting d^TZ from d~TZ and if 
there is no superfluous edge in /. This notion of interface corresponds to the geometrical notion 
of interface if, to the edges of / we associate their dual d— \ dimensional facets. 

According to the max-flow min-cut Theorem [7], the maximum flow from d~TZ to d^TZ by the 
edges of capacities Je is also the flow through the interface of minimal capacity. We use this 
characterization for our deflnition. Given a rectangular parallelepiped TZ — TZa:^L,HiS,n) as in 
()1.15p we call I{TZ) the set of interfaces for TZ and define the maximal flow in TZ, for a realization 
J of the media, as 

(1.28) Mi = T:^ inf Vje. 

eel 

This quantity has the same properties as surface tension since it is as well sub-additive. In 
particular, the maximal flow in TZo,N,SN{S,n) converges in P-probability, upper deviations occur 
at volume order and lower deviations occur at surface order [411 US] ■ We will make use of the 
following results: 
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Theorem 1.10. There exists fJ-{n) G [0, +cxd), the maximal flow for the distribution P in the 
direction n G S"^^^ , such that, for any 5 > and S G §ri, 

f^Tlo.N.sN{n,S) j^:^ Kn) in P-probability. 

It is positive if and only if¥(Je > 0) > Pc{d). Furthermore, 

(1.29) J"""||n||i s^M(n) 
and the inequality is strict when P(Je > J™'") > Pc{d). 

The convergence of the maximal flow is a consequence of the sub-additivity. It is shown in [46j 
that the maximal flow is when P(Je > 0) ^ Pc{d), while its positivity was established in [T3], 
under the conjecture that the critical threshold for percolation and slab percolation do coincide, 
proved later on in [26 . The inequality (|1.29p is easily obtained from the remark that minimal 
interfaces have cardinal of order A^'^"^||n||i. When P(Je > J™'") > Pc{d), for small e > there 
is a percolating net of edges with values Jg ^ J""" + e [26l [38] , which is responsible for the strict 
inequality. See also Proposition 4.1 in 39J. 

It turns out that the maximal flow determines the asymptotics of the quenched value of surface 
tension at low temperatures. Precisely, we show that; 

Proposition 1.11. LetF be a product measure on [0, l]'^ such thatF{Je > 0) = 1. Then, uniformly 
over n £ S"^'-^ , 

(1.30) hm ^ - /.(n). 

Clearly, p.30p also holds in the case P(Je > 0) ^ Pc{d) since T0{n) ^ /3/i('n.) — (cf. Lemma 
13.11 and Theorem I l.lOp . When P(Je > 0) > Pc{d) a renormalization argument allows us to prove 
that: 

Proposition 1.12. Assume that P(Je > 0) > Pc{d). Then, 

(1.31) liniinf^^>0, 

/3^+oo P 

uniformly over n G S'^~^ . 

On the other hand, the surface tension under the averaged Gibbs measure is asymptotically 
determined by J""": 

Proposition 1.13. For all product measure P on [0, l]'* and all A > 0, uniformly over n G S^^^ , 

(1.32) hrn l£^ ^ AJ-'"||n||i. 

In the case J""" = and P(Je > 0) = 1, an equivalent to Ta{n) is given in Proposition 13.21 
These asymptotics have consequences on the shape of the crystals under both the quenched and 
the averaged Gibbs measure (see Proposition 13.31 below) . They also immediately imply that the 
inequality r^ (n) ^ Ar? (n) is strict at low temperatures in a number of cases: 

Corollary 1.14. Assume that P(Je > J™™) > Pc{d). Then, for any A > there is (3^ < oo such 
that 

(1.33) T^(n) < AT«(n), Vn G S'''-\V/3 > /3^ 

In particular if J""" = and if there is still a phase transition (i.e. P(Je > 0) > pc{d)), then 
the inequality is always strict at low temperatures. 

One consequence of (|1.33|) is the strict inequality Tj^™{n) < T%{n) under the same assumptions 
than in the Corollary, as ArS""(n) ^ T^in) (Proposition I2.3p . 

Let us conclude on a comparison with the directed polymer model in 1 + 1 dimensions: for the 
latter model, it was proved recently T7J that the Lyapunov exponent is positive at all /3 ^ 0, which 
corresponds in our settings to the strict inequality TS(n) = '^b^^ < '''ai''^)- 
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1.5. Phase coexistence. We describe finally the phenomenon of phase coexistence in the dilute 
Ising model. Phase coexistence occurs when both the plus and the minus phase are present at the 
same time and occupy (distinct) regions of the domain. This phenomenon does not occur naturally 
in the Ising model. One way of obtaining phase coexistence is by conditioning the measure /x^ on 
the event that the overall magnetization 

(1-34) "^A = -r^Xl^^ 

is smaller than m < mp. Under this conditional measure, we will show that the two phases do 
coexist and that the minus phases occupies a fraction of the volume v such that (1 — 2v)mp = m. 
Furthermore, the shape U of the region containing the minus phase is deterministic : if r is the 
surface tension of the model, the observed shape minimizes the surface energy 

T{U) = I T[n)ds 

JdU 

under the volume constraint Vol(C/) ^ v, and this implies that [/ is a translated of v^^'^W where 
W is the renormalized Wulff crystal associated to t: 

(1.35) W = \{x(^M.'^ ■.x-n^T{n),\fne S"^-^} 

where A > is chosen such that Vol(yV) = 1. 

Before we state our results, let us recall that A// stands for the at-most-countable set of /? at 
which lower large deviations for surface tension are possible at less than surface order (Corollary 
and (3c is the slab percolation threshold (|1.14p . Another important notation is 



(1.36) 7V= 1/5 ^ : lim E$f'^ ^ lim E^f'") 

the set of /3 such that infinite volume averaged FK measures are not unique. J\f is at most countable 
(Theorem 2.3 in 44J). 

We denote by W (resp. W^) the Wulff crystal associated with the surface tension t'^ (resp. 
T'*') as in (jl.35p . and v — a'^ the fraction of the volume occupied by the minus phase. The Wulff 
crystal aW of volume v fits into the unit box [0, 1]'' only if Q;diamoo(W) ^ 1, where 

diamoo(^) = sup \\x~y\\^, AdW^. 

x,y^A 

Our first theorem concerns the cost of the lower large deviations for the magnetization. In the 
sequel, A^ = {1, . . . , iV}'*. 

Theorem 1.15. Assume (3 > J3c with /? ^ TV. Then, for all ^ a < l/diamoo(W'^), 

(1.37) ^-rlog^'^'+ ( ^^^^ s; l-2a'' I — >-J^''(aW«) mf -probability. 

Then we describe the geometry of the two phases. We consider a mesoscopic scale X G N* and 
define the magnetization profile A^i<- as 

(1.38) ^-^ [«'l]' - t-^'l] 
where 

(1.39) i{x)^[ ^ ,..., ^ ) and ^, ^ Ki + {I, . . . ,KY 



'W l^zeANnAi^^i ^z 



Hence, unless x is too close to the border of [0,1]'*, M.k{x) is the magnetization in a block of 
side-length K that contains Nx. Theorem 5.7 in [33] provides a strong stochastic control on Ali^: 
when (3 > [3c- In particular, when K is large enough, at every x the probability that M.k{x) is 
close to either mp or —mp is close to one under the averaged measure E/i^' . Hence Mx/'mp 



'Nxi' 
K 


,..., 


'Nxd 
K 
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describes the geometry of the phases in the Ising model: it is close to one on the plus phase region, 
close to minus one on the minus phase region. 

We need a few more notations. To U cW^ Borel measurable, we associate the profile 



(f.40) 



Xc/ : a; e K i-^ 



1 if x I 
— f else 



U 



We also consider the set of vectors 



and denote by ||.||li the norm of the L^-space L^ ([0, 1]''; 
z such that the translate z + U fits into [0, f]*^: 

(f.41) T{U) = {zeR'^ : z + U (l[0,lf}. 

Our second theorem describes the geometrical structure of the two phases when they coexist: 

Theorem 1.16. Assume that (3 > J3c and (3 ^ M . For all ^ a < l/diamooCyV) and e > 0, for 
any K large enough one has 

€ 1 - 20" = 1 



(1.42) lim ^;(^+ ( inf 


Mk 

mp 


^ Xz+aW 




mp 


in ^-probability (T-a.s. when /3 ^ A// 


)■ 









Note that, although we state our theorems for the Ising model, they could easily be adapted to 
the Potts model with random interactions, or to random-cluster models, as the two fundamental 
tools for the study of phase coexistence, the coarse graining [44] and the study of surface tension, 
were developed in the more general setting of the random-cluster model (g ^ 1) with random 
couplings. 

The fact that we consider K large but finite is a slight improvement with respect to former 
works. In general, one can take any K = K^ such that 1 ^ Kn ^ N because on the one hand, 
Mkn is close to the local mean oi Mk as Kjy » 1, and this local mean is close to Xz+aW because 
the iiTAT-blocks intersecting Nd{z + aW) contribute to a negligible volume as Kn <C N. 

Let us conclude this paragraph on a first consequence of the presence of random couplings : 
the limit shape of the droplet at low temperatures is smoother. In the case of the pure Ising 
model, the Wulff crystal converges to the unit hypercube [±1/2]'' as the temperature goes to zero. 
Here, W"? converges to the Wulff crystal associated with the maximal flow ^ when P( Je > 0) = 1 
(Proposition [131) . Little is known on the crystal VV associated to the maximal flow ^. Yet, as 
discussed in Section [3^ an argument by Durrett and Liggett [H] shows that W^ is not a square 
when, for instance, d = 2, 

P {Je = 1/2) = p and P ( Je = 1) = 1 - p 

with pc < p < I, where Pc is the critical threshold for oriented bond percolation. 

1.6. Phase coexistence under averaged Gibbs measures. Now we consider the issue of phase 
coexistence under averaged Gibbs measures, that is, when phase coexistence is imposed on both the 
spin configuration and the random couplings. Before we go further, we would like to remark that 
averaged Gibbs measures do not have the physical meaning of the quenched measure: in quenched 
ferromagnets, the disorder is frozen and thus cannot be influenced by the spin configuration itself. 
However, the analysis presented here gives an insight on the phenomenon of localization which can 
occur in models with media randomness. 

First we remark that the cost for phase coexistence is here determined by the surface tension 

Theorem 1.17. For all X > and ^ a < l/diamoo(yV^); 

A 



(1.43) 



ATd 



— logE 



^^" * mp 



N- 



-T^{aW^). 
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The inequality r'^ < At' at low temperatures (Corollary 1 1.14|) implies that 

in other words the cost for phase coexistence is strictly smaller under the averaged Gibbs measure 
than under the quenched Gibbs measure. One can go further and analyze the cost for reducing 
the cost for phase coexistence under averaged Gibbs measures: the functional 

j7(/) = sup{.F^(W^)-A/}e[0,oo], /eR 

A>0 

is the rate function for lower deviations of the cost for phase coexistence. If W™'" and jF™™ stand 
respectively for the Wulff crystal and the surface energy associated to r™'", then J' is infinite on 
the left of jrmin(-y^;min)^ g^^j^g ^^ ^^le right of jrmin(>vmin) ^^^ ^cro On the right of J^«(yy«), and: 

Corollary 1.18. For any f ^ jrmin(yvmin) and a ^ small enough, 



lim 



1 



N N"^- 



-logl 



/ 1 



ViV''- 



llog/^A 



J,+ 






^ 1 - 2aM ^ -a'^-'f = 



''Jif)- 



Upper deviations for the cost of phase coexistence, on the other hand, happen at volume order 
(cf . the proof of Proposition I4.14p . 

The shape of crystals under averaged Gibbs measures is as well determined by the surface tension 

T^(n): 

Theorem 1.19. For any ^ a < l/diamoo(yV'^^^) and e > 0, for any K large enough one has 



(1.44) 



lim 



^Ma« 



inf 



M 



K 



mp 



Xz+a'W^ = '^ 



^e 



Li 



nip 



€ 1 - 20"* = 1 



This result extends in fact to all A > 0, at the price however of heavier notations, because 

A, 



'E[{fjLj^'^{.))^] is not a measure when A ^ 1 



Theorem 1.20. For any A > 0, any ^ a < l/diamoo(W^) and e > 0, for any K large enough 
one has 

A" 



E 



lim ■ 



J,+ 



inf 



z&T{aW^ 



Mi 



Xz+aW>- 



< e and 



^ 1 - 2a' 



E 






^ 1 - 2a' 



= 1. 



We conclude the summary of our results with a description of a phenomenon of localization. 
First, let us characterize the typical value of the surface tension r^ under the averaged measure, 
conditioned to phase coexistence. For any /3 ^ 0, A > and n G S'^^^, this value stands between 

(1.45) T^'-{n) = inf {t ^0 : /„(T) + AT = T^(n)} 

(1.46) and f^'+(n) = sup {r ^ : /„(t) + At = T^(n)} . 

The equality T^^~{n) = T^^'^{n) holds whenever there is at most one r at which the slope of /„ 
equals A, that is, for all but at most countably many values of A > 0. Note also that the strict 
inequality T^{n) < Xt'^{ti) implies f'^'+(n) < T'^{n). We also consider a similar quantity for the 
quenched value of surface tension: 

(1.47) f«(n)=inf{T:/„(T)=0} 

which coincides with T'^{n) for all n E S"^^^, for all but at most countably many /3 J^ 0, see 
Corollary 

Our last Theorem describes the value of the surface tension r'^ conditionally on the position of 
the crystal: we prove that the typical value of surface tension is r* outside the boundary of the 
crystal, while on the boundary it is reduced to f'*'. When T'''{n) < XT'^{n) for some n G S"^"^, the 
location of the Wulff crystal under averaged Gibbs measures is thus determined by the realization 
of the media: the boundary of the crystal coincides with the place where surface tension is reduced. 
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Theorem 1.21. Let A > 0, ^ a < l/diamoo(>V^) and z e T{aW^). Consider h,6,^ > and a 
parallelepiped rectangle TZ — TZxji.5h{n,S) C (0, 1)'' as in hi. 15]) . Call TZ^ — NTZ + zn{TZ) where 
zn{TZ) € (—1/2, 1/2]'' is chosen such that the center ofTZ^ belongs to Z''. For e > small enough 
and K large enough, 



E 


'{i^i: 


T.^jv e A and 


Mk ^ 

ran ^ 


z+aW^ 


.-))1 


\f^oo 


E 


(^t^ ( 


Mk 


Xz+aW^ 


.-^))1 





and A = [f«(n) - 7, riin) + 7] 



n IS 



the outer local normal to z + aW^ at x, h is small enough and 



when 

i. nnz + adw^ = a 

ii. or X d^ z + adW^ 

A= [f^'"(n)-7,f^'+(n) + 7] 
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PhD Thesis at Universite Paris Diderot [IS]. I am also grateful to Marie Theret and Raphael 
Rossignol for useful and pleasant discussion about maximal flows and concentration. 



2. Surface tension 

As announced in the former Section, surface tension is a fundamental tool for understanding 
the mechanism of phase coexistence. It quantifies the free energy per surface unit of an interface 
separating the plus and minus phases in the dilute Ising model. In this Section, we prove the 
convergence of surface tension in dilute models and study its large deviations. 



2.1. Sub-additivity and convergence. In many aspects the surface tension for the dilute Ising 
model is similar to the one of the Ising model with deterministic couplings. It has the crucial 
property of being sub-additive, as in the uniform case |36j : this is shown in Theorem 12.11 below. 
We present here the proof of Proposition 11.21 Theorem 12.11 and finally Theorem 11.31 We also 
explain why surface tension is positive under the assumption that f3 > f3c (Proposition 1 1 . 5| ) . 
Proof (Proposition 1 1.2p . The surface tension x^ is a non-decreasing function of J and /3 because 
D-ji is a decreasing event while the measure $.^ stochastically increases with p = 1 — exp(— /3Je). 
Now we consider H' ^ H and call TZ = TZx.L.H{S,n) and TZ' = TZx,L,H'{S,n). In view of the 



^J,7r 



J.w 



DLR equation and of the monotonicity of ^^ along tt, the measure ^^, restricted to E{TZ) is 



\J,W 



stochastically smaller than ^^ . On the other hand, it is clear that V-jz C T>tz', and because V-jz 
is a decreasing event we conclude that 



""" ^ Tj^ ^ T^'^'^ is a consequence 

max 



which shows that tj^ is a non-increasing function of H. 

It is clear from the definition that t.^'" ^ 0. The inequality t^ 
of the monotony in J. We conclude with the upper bound on t^'^^. Because of the monotony in H 
we can take H — 2\fd (which ensures that disconnection is still possible). We have: t^'^ ^ ^v?^ 
where TZI — TZ^ ^ 2\/d('^' '"')■ ^* ^^ enough to close all the edges of TZ' to realize the disconnection in 
TZ' . The DLR equation, combined with the monotonicity of <& 
TT yields: 



'K' 



^-7^i°g n 



$ 



{e} 



({w, = 0}) = /3J 



r 1 ^ along the boundary condition 

\E{TZ')\ 



Ld- 



eeE{TZ') 
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Finally, \E{TZ')\ is not larger than 2d times the cardinal of 7?.', which is itself not larger than the 



volume of V^ = y 



xeTt' 



[0,1]'') cn„j^^2Vd,3Vdi^^'^)- Consequently, 



'K' 



^ /3 J"^'' X 2d X 



(i + 2Vd) ' 



eVd 



j^d-i 



< /SJ""'*'^ X 6 X 2^*^3/2, 



D 



Now we address the issue of sub-additivity. It is a fundamental tool not only for proving the 
convergence of surface tension, but also for establishing the large deviations principles in the next 
Section. 



Theorem 2.1. Consider n e S'^^^ , S,S' C §„ and H,l ^ 2Vd, L > iVdl. 
TZq ^ H+\/~d/2^^^''^)- There is a collection {TZi)i£c of rectangular parallelepipeds TZi ~ TZ^ 



Let n = 



that are disjoint subsets of TZ, centered at Zi £ Z 

I r 



(2.1) 

such that, for any J : E{TZ) 

(2.2) 



i-Cd 



[0,1]; 



<^ 



L I 



< 



' ' lec 



with 

I 
L 



(3cd 



d-l 



\C\^l 



where Cd < oo is a constant that depends on d only. 

Let us make a few comments on this Theorem. First, a key feature of the sub-additivity as 

formulated in Theorem 1 2. II is the independence of the r^. under P since the TZi are disjoint. Note 

that as well, the x^. have the same law as the TZi are all centered at lattice points. Three error 

terms appear in Theorem 12.11 Their origins are as follows (see also Figure |31): 

i. the term j3cd/l stands for the cost of disconnection in the middle section of TZ between adjacent 

ii. the term (3cdl/ L represents the cost of disconnection in the area not covered by the TZi 
iii. and the increase of H by V"/2 for TZ with respect to the TZi is a consequence of the requirement 

that the TZi be all centered at lattice points. 
The last error term could be avoided for rational directions n e S"*"^, yet (as the two others) it 
will soon disappear when we take the limit 7f — > oo. 




Figure 3. The rectangular parallelepiped TZ and the collection {TZi)i^c in Theorem 12. II 
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The reader will notice that the use of the FK representation permits to give a relatively short 
proof of Theorem 12. II 

Proof (Theorem 12. ip . We begin with the definition of Zi and C. We call {e'^)k=i...d-i the edges 
of 5' and e^ — n, so that {e'i^)k=i...d is an orthonormal basis of M"^. For all i = {ik)k=i...d-i G Z"^'^ 
we define Zi as the unique point of Z such that 



z. 

k=\ 

and call 

7d-l 



2' 2 

c = {i e if-^ : n., c 7^} 



letting 

"^ = '^o,L,H+Vd/2('^'"') s-i^d TZi=nz^^LH{S',n). 
We proceed with the proof of (j2.1|) first. We call 7i„ the hyperplane of M'' orthogonal to n that 
contains and remark that the orthogonal projections of Zi + IS' (for all i g C) on Tin are disjoint 
and all included in LS. Hence their total surface |C|Z''~^ does not exceed the surface of LS, namely 
L'^~^, and the upper bound in (|2.ip follows. Reusing the previous notations we call 

d-l 

z^ = (i + Vrf) ^ ike'k, Vz e Z'^-^ 

k=l 

so that z,' G Tiri- We consider then 

C = Ue IJ^-^ : z', + (l + Vd\ S' C LS\ . 

In view of the inequality d{zi, zQ ^ y/d/2 it follows that Zi + IS' C TZ, for all i g C, hence C C C. 
On the other hand, for any i e Z'^^'^ such that z,^ + (? + \/rf)5' n (L - 2\/d(Z + Vd))S 7^ we have 
i G C, hence 



and 



pd-i 



l + Vd 



:^ 




^ \-{d-\){^^2^d^- 



which yields the lower bound for ()2.ip . We pass now to the proof of (12. 2p and call 



2 



where d{e,Ti.n) stands for the shortest distance between one extremity of e and Tin- The inclusion 



\iec / 



holds: consider to that belongs to the left-hand side and let c an w-open path issued from d^TZ. 
Every times c enters some TZi by the upper boundary 9+7?,, it also exits by the same upper boundary 
since w S V- . As c cannot use the edges of £ it is not able to cross the middle hyperplane Tin 

elsewhere than in the TZi, and in particular it cannot reach d^TZ. Since the T)^ as well as the 
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{tJe — 0} are decreasing events, the DLR equations and the nionotonicity along the boundary 
condition for ^'^ imply that 



(2.3) ^ n'^^r(^K.)xexp(-/?ifi) 






as $,^7({a;e = 0} = 1 — Pe = exp (— /3Je) ^ exp(— /?). We proceed then with an estimate over the 
cardinality of £: we call F = {a; G Z'' : 3j/, {x,?/} G £} the set of extremities of some e G £ and 
remark that \£\ s: d Vol [V) where y = U:EeF ^^ + [0, 1]''. We have 

y ^'^^.L+2Vd,:^Vd/2iS,n) while ^ n 7^,^,_2^_^(5^n) = 0, V* G C, 
hence 

in view of the lower bound in ()2.ip . Taking logarithms in (j2.3p and dividing by —L'^^^ we obtain 
the inequality 



ri^fiV^Vr, 






r4+'^./^(^ + j) 



and (|2.2p follows from the upper bound in (|2.ip . 

We conclude with a word on the structure of the sequence {rj^JieC- The 7?.i are disjoint by 
construction, hence so are the edge sets E{iti), hence the tj^. are independent. They are identically 
distributed as the TZi are all centered at lattice points, P being translation invariant as a product 
measure. D 

Now we establish the convergence for surface tension and prove Theorem II .31 The proof of this 
Theorem is based on the sub-additivity of surface tension. We do not apply directly Kingman's 
sub-additive Theorem [34j as we want to show that t"? does not depend on S, nor on S. 
Proof Taking the expectation E in the sub-additivity inequality (|2.2p we get 



Applying linisup^^g^, then liniinf;^oo and taking the decreasing limit in H we obtain 
lim limsupEr^ (•<;^> ^ lini liminfEr^ (•<;' .^-i 

which proves that 

(2.4) T'^(n) ^ lim liminfEr.^ r<;^\— hm limsupEx^^ i<;^\ 

exists and does not depend on 5 G §„. 

We prove now the convergence t^n -^ T'^{n) in P-probability, where Ti^ = TZq, n,sn{S, n). The 
sub-additivity (12. 2p yields: for any 5 > Q and N large enough, 

^K« ^ ''k>,N,„+Vl/2(S:r^) =^ j^ Zl ^K., ,z,.« + P^d (jj + J 

Taking limsup^^^^ and applying the strong law of large numbers give: 



limsupr^iv < Er^o^^(5,„) + -p IP-a-S. 

and after liminfi_>oo and lim//^oo we see that, for all 5 G Sji and (5 > 0, 

(2.5) lim sup T^w s;T«(n) P-a.s. 
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On the other hand, the sub-additivity (12. 2p is also responsible for the convergence of Et^n : remark 
that 

Et-^ ^ (5 ») ^ ^^n« + /^Cd ( — + — 

hence limsup^^g^^ followed by lim inf Ar_>oo give: 

(2.6) T''in) =^ liminfET.^iv. 

Together with ()2.5p and (12. 6p . the boundedness of t.^„ ensures the convergence in probability. D 

Let us sketch now a proof of Proposition 11.51 namely that the quenched surface tension r"? (n) 
is positive for any (3 > (3c'- thanks to the renormalization argument of j44j . one can compare the 
surface tension r" — t^^^ under the averaged Gibbs measure to the surface tension of high density 
site percolation, which is positive. The claim follows as r* ^ t° by Jensen's inequality. 

2.2. Upper large deviations. Due to the presence of the random couplings, surface tension can 
fluctuate around its typical value. The sub-additivity permits to study the order of the cost of 
large deviations. First, we examine upper deviations and prove Theorem 11.41 The proof is based 
on the following argument: we split TZo^N,5N{S,n) into cN rectangular parallelepipeds TZi with 
finite height H. In order to increase r^ /c \ one has to increase surface tension in each 7?,,, 

but the cost of increasing one tj^. is already of surface order by sub-additivity. 
Proof f Theorem II .4^ . As a first step towards the proof we estimate the cost for upper deviations 
of surface tension in a rectangular parallelepiped of fixed height, using the sub-additivity of t"^ . 
From the definition of T''{n) at (I2.4p it follows that for any H large enough. 



limsupET7j^^^^(5 „) < T^in) + 



Given such an iJ we fix Z large enough such that Et^ ,^ ^ ^ r'^iri) + e/3 and CdP/l ^ e/4, 
where Cd refers to the constant in the sub-additivity equation. With the notations of Theorem l2.1l 
we have: 

and the r^ (s n) ^'^^ i-i-d. variables of mean not larger than T'^{n) + e/3. Hence, Cramer's 
Theorem tells that 




£ 



,l,H 



(5,r^)^^'(") + ^ ^exp(-c|C|) 



for some c > 0. Reporting in (|2.7p proves that for any e > 0, for any H large enough: 



1 



(2.8) linasup--^logPfr^^,^^(5,,) ^ r^n) + e < 



L^oo J^ 



- that is, the cost for increasing t^ .^ . is of surface order. We fix such an H and decompose 
now the rectangular parallelepiped TZ = TIq, n.sn{S, n) in the direction n. Precisely, we let 

ii = 2Lff+— Jm, Vi e Z and Ht = U^^j^jj^^/^{S,n). 

We call G the set of i G Z such that TZi C 7?. and consider, for all i E Q, Xi the point of Z'' such 
that Xiex^ + [-1/2, 1/2)^^ and let 
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The rectangular parallelepipeds TZi are disjoint subsets of TZ = TZo,N,SNiS,'n), all centered at 
lattice points. Furthermore, if we call £iat the set of edges in E{TZ) with one extremity at distance 
at most V" from the lateral boundary of 7?., we have: 

w S M V-ji- and uje = 0,Ve G £iat => ^ G 2?k- 

Hence the DLR equation yields: 

$^"'(I?k) ^ max$^^"K=0,Vee£iat and u; G Pk J 
^ g-'^l^'-'l X max^^"" (w G PkJ • 

As |£iat| =^ CdSN'''^^ we conclude finally to the inequality 

(2.9) T-^ ^ CdSf3 + mmT^.. 

ieG ' 

Inequality (j2.9p states that in order to increase significantly r^^, one must increase each tj^.. Yet, 
the cost for increasing one of the x^ . is of surface order (|2.8p , and the x^ . are independent variables. 
Hence for any S > such that Cd6/3 < e, 

hmsup -1^ logP (r^^ ^„(5^„) ^ r«(n) + 2e) < 0. 

As T^ /ON decreases with (5, the claim follows for arbitrary S > 0. D 

2.3. Lower large deviations. Contrary to upper deviations, lower large deviations occur at 
surface order. Here we consider the rate function /„ for lower large deviations. The fact that 
deviations occur at the same order as the disconnecting event defining surface tension is responsible 
for the distinct behavior of surface tension under quenched and averaged measures. Explicit bounds 
on the rate function /„ will be derived in Sections 12.51 and 12.61 

Proof We begin with the definition of the rate function I-ji in a rectangular parallelepiped 
7?. = 7?.o,L,H('5,n) as the surface cost for reducing tj^ to t: 

lTz{T)=-j^l0g¥{Tl^T). 

According to Proposition II. 21 I-ji^^ ^ His.n) (''") is a non-increasing function of r and H . Hence the 
limit 

(2.10) /(5^„) (r) == lim inf limsup/^^ ^ « (5.n) (t + e) G [0, oo] 

exists - we introduce the parameter e > in order to compensate for the error terms in (12. 2p . It is 
clearly a non-increasing function of t. We prove now that it is also convex in r and that it does not 
depend on 5 G Sn- let S' G §„, e > and a G [0, 1]. Using the notations 7?. = 7?.q ^ H+y/d/2(^^ '^)' 
TZi — TZzi.i.H{S' ,n) and C of the sub-additivity Theorem (Theorem 12. ip . we have 

if C^ U C^ is a partition of C such that 

(2.11) rl^l^l^' ij^^g 
The probability for realizing condition (j2.1ip equals 

e^p{^\CY-'In„,MS',r^) {r' + e) - \CY-'lno,MS'.r.) {r' + e)) 
and letting |C^|/|C| -^ a and i -^ oo we see that 

(2 12) limsupi ^T^o.L.H+^.M^^^) ("^' + (1 - ")^' + 2e + (3cd/l) < 

aIno.i.H{S'..r,) {t^ + e) + (1 - a)/^^ , ,^(5-,„) (r^ + e) . 
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Taking the superior limit in I, then the hmit in H, then e ^ 0+ we obtain 

I{S,n,) (ar^ + (1 - a)T^) ^ al^g-.r,) (r^) + (1 - a)I(s',r,) (r^) 

which proves both the independence of I(s,n) with respect to S (take a — 1) and the convexity 
along r. We let now /„ = /(5.„) and postpone the proof of (jl.lSp for a while. The continuity of 
In on the interior of the domain of finiteness of /,,, is a consequence of its convexity. Hence we 
examine the domain of finiteness of In- Let first r < T"^"^{n). If e > is small enough, the event 
T^ (^ .^> ^ r + £ < T'^™{n) has a probability zero and consequently, In{T) = +oo. The second 
easy regime is r ^ T'^{n): from Proposition II .41 we infer that limi^oo ^{t^ rs n) ^ ^ + ^) ^ -'^ 
provided that H is large enough and this implies In{T) = 0. If at last r > T™^^^{n), there is H 
such that 



limsupr7^""^(5^„) <T. 



L 

We will prove that, for 5 > small enough we still have: 
(2.13) limsupr.£""°+'',<, , < r. 

If we let TZ = TZo^l,h{S, n) and differentiate along S, we obtain 



dS ^ dJe 

eeB(K) 



,7=J"""+5 



yet, (|2.22p and Proposition 12.61 indicate that for any J ^ J ^ 



d-i a^j 
n 



U'-^ dr. 



P dJe 



s$ 1. 



As a consequence, r^ is a C(j/37J-Lipschitz function of S. The same is true for lim sup^ r^ (s n)' 

thus (|2.13p holds true for (5 > small enough. Now we write, for any L large enough: 



WL,H{S,r.){T) = — _-^logP(r^^^^(5^„) sc:r 



^ --^logP(./eS=; J'"'" + 5, yeeE(TZo,L,H{S,n: 

< CdH X (- l0gP( Je e [J""", J""" + 5])) 

which is finite thanks to the definition of J""". This ends the proof that In{T) < oo, for any 

We address at last the convergence (|1.18p . The inequahty lKo,N.5N{s,n) (t) ^ /7?,o,N,ff(5,Ti) (t) 
when iV(5 ^ iJ yields an upper bound on the superior limit: 

lim sup /-K^ „ sNiS,n) (t) ^ inf limsup/K,, ^ ^(s.n) (t) s$ /^(t") = /„(r) 

for all T > T™™(n), thanks to the continuity of /„. For the lower bound we use the sub-additivity 
of surface tension. Applying (|2.12p with a = 1, I = N, H ~ SN yields: for any e > and N large 
enough, 

lim^sup/7j^^^^^^^^(5,,^) {r + 3e) < lTZo.^,sN{S.r,} {t + e) 
and replacing t + £ with t, we obtain after the limits iV — > oo and e ^ 0+ the lower bound 



/„(r) ^ liniinf /k„ „_,„(5,„) (r) , Vr G 



n 
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2.4. Surface tension under averaged Gibbs measures. The rate function /„ can be analyzed 
througli a dual quantity: the surface tension under the averaged Gibbs measure defined at p.24p . 
The duality of Fenchel-Legendre transforms for convex functions (Lemma 4.5.8 in [18j ) implies 
that A 1-^ T^{n) is concave and that 



(2.14) 



/„(t) =sup{T^(n)-Ar}. 

A>0 



As we said at (|1.25p . T^{n) can be interpreted as the surface tension under an average of '^■j^ ■ 



Indeed, if we let 



(2.15) 



'n 



Ld- 



-logE 



'^'tT i-Dn) 



= --^^logE(exp(-AL'^-V^)), 

for any rectangular parallelepiped TZ of side-length L as in (jl.lSp . then Varadhan's Lemma yields: 

Proposition 2.2. For any A > and n G S^^^ , for any sequence of rectangular parallelepipeds 
Ti.^ — TZo^n,sn{S, n) with 6 > Q and S E Sn, the quantity t^n converges to T'^{n): 



(2.16) 



lim T, 



N 



K" 



r\n). 



Thus, the limit does not depend on 5 > Q nor on S E Sn- 



We defined at (|1.47p the value f'(n) of the surface tension at which /,i(t) becomes zero. Below 
are some immediate consequences of the definition of t'^ (n) at (|1.24p together with (|2.16p , which 
allow to sketch the graph of A i-^ r'*' (n) on Figure UJ 

Proposition 2.3. The following inequalities hold: 

(2.17) AT™'^(n) =^ T^(n) =^ Af'(n), Vne5'^-\A>0 

while: 



(2.18) 



A A^0+ 



f'(n) and 



Hn) 



A A^+c 



T'"'"(n), Vn e S'^^K 



Hence, T^{n) is positive if and only ifT'^{n) > 0. Furthermore: 



(2.19) 



T^{n) — > lim /„(t) e [0,oo]. 



/n(()) 



» Af'(n) 





Figure 4. The graph of A i-^ T^{n) in the case of dilution (r™'" = and /„(0) < 
oo, left) and distributions with t™'" > (right). 

Another important yet classical fact is the convexity of surface tension [33]. The proposition 
below is a consequence of the weak triangle inequality for tj^ (see [3S] or [10] for the uniform case, 
or Appendix 2.5.2 in [i5]). 

Proposition 2.4. Let /' be the homogeneous extension of t"^ to W^, namely: 

||a;||T9(a;/||a;|l) ifxER''\{0} 
ifx^O, 



n^) 



and let f^ (resp. f^) he the homogeneous extension of t^ (resp. f^) to 
are convex and r^, t'^ and f are continuous on S'^^^ . 



Then, /«, /^ and /« 
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2.5. Concentration at low temperatures. In this Section and tlie next one we establish 
respectively Theorems 11.71 and 11.81 In both cases we use concentration of measure theory, which 
is a very efficient tool for analyzing the fluctuations of product measures. In the case of polymers 
or even spin glasses it yields relevant bounds on the probabilities of deviations, see [35j for a 
review. Concerning the Ising (or random-cluster) model with random couplings, its application 
to the deviations of surface tension requires a control over the surface of the interface, and this is 
the point where the proofs of Theorems 11.71 and 11.81 differ: at low temperatures one can control 
rather easily the length of the interface, while under the only assumptions of Theorem ll.Sl the same 
control is not immediate. 

The surface tension T'^{n) under averaged Gibbs measure plays an important role here, as well 
as the modified measure E\ defined at ()2.25p below. We will obtain lower bounds on T''^{n), which 
correspond to lower bounds on Inir) by (j2.14|) . 

Rather than making the assumption that the product measure P satisfies a logarithmic Sobolev 
inequality as in 45|j, we use general bounds on product measure (Corollary 5.8 in [35] ). The 
author thanks Raphael Rossignol for pointing out this improvement. The proof of Theorem II. 71 is 
made of four steps, the first three being common with the proof of Theorem II .81 

The first step consists in relating the derivative of the surface tension r^ (n) in a rectangular 
parallelepiped TZ as in (|1.15p . with a basis of side- length L, to the entropy of the positive function 
exp(/A) where 

h = -AL'^-V^. 
We recall that the entropy of a positive measurable function / with E(/log(l -|- /)) < oo is 
(2.20) Entp(/) =E(/log/)-E(/)logE(/). 

With these notations, it is immediate that: 



Lemma 2.5. For any A > 0, 




(2 21) -' (^^) = 
^ ' dX\\ ) 


1 Entp(exp(/A)) 
A2L'i-i E(exp(/A)) 


As a second step we study the quantity 




(2.22) ai = 


L-^-i dr:^ 



(i dJe 

The proposition below provides an interpretation of a^ as the probability that the disconnecting 
interface due to the event I?k passes through the edge e. We prove also, and this is crucial for our 
construction, that the actual value of Je does not influence too much that of a^ : 

Proposition 2.6. For any e, a'l is a C°° function of the Je' ■ For any J e [0, 1]^(^ \ one has 

(2.23) a,^ = -($^"(c^e)-$K"(^el^K)) tf Je>0 

together with the following inequalities: 

(2.24) ^ af ^ 1 and supaf < e'^ inf af . 

The controls (|2.24p . together with Corollary 5.8 in [3S], permit to establish the third step. Given 
a rectangular parallelepiped TZ as in ()1.15p and A ^ 0, we introduce the probability measure Pa 
that to any bounded measurable h : J ^^ h{J) G I^ gives expectation 

/ exD i—\L'^~^'^'^ 

(2.25) Ea(/i(J))=eU(J) ^^ 



Eexp(-AL''-iT7^ 




Usual measures such as dilution P(Je G {0, 1}) = 1, or Je with positive density on [0, 1] do satisfy a logarithmic 
Sobolev inequality, cf. |35j or Theorems 4.2, 6.6 and Section 6.3 in |12) . 
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Proposition 2.7. Denote mp = E(Je). For any A ^ 0, we have both 



(2.26) 



Entp(exp(/A)) 
E(exp(A)) 






2p/3(l+A) 



Ea E 



m;3 A 3/3 



The second majoration leads to Theorem II .81 while the first one yields Theorem 1 1 . 71 after a last 
step: using Peierls' argument we show that, in the Ising model {q = 2) with couplings Je ^ e > 0, 
the length of the interface is of order N'^~^. 

Proposition 2.8. Let q — 2 and e > 0. There exists c^ < oo such that, for (3 large enough, for 
TZ = TZ^ = TZ().N.SN{S,n) with S e (0, 1), for any realization J of the random couplings such that 
Je ^ s and N large enough, 

(2.27) Y. «e «; ^TV^-i. 

We give now the proofs of all the propositions, followed by that of Theorem 1 1.71 
Proof (Proposition 12. 6p . The fact that a^ is a C°° function of Jg' is a consequence of the same 
property for tj^, the quantity $^"'(Pk) being always positive. We introduce next a few notations: 
we let 

Pe 



(2.28) 



wni^) = n 



for any w e ^e(tz) 
we have 



eeE{n) 
''E(n) 



1-pe 



qCZm^-) and Zi{A)^Y.<^^) 



ujEA 



and A C r2p,T3^, see (|1.9p for the definition of C"" - (cj). For aU J with Je > 0, 



EilZ)'' 



d\ogw^{uj) ^^^ 

dJe Pe 



and as a consequence, for all J with Je > 0, 



P. 



(3 dJe Zj^{n^^^s^) 



.J,w , 



Under this formulation, the FKG inequality and the bound $^™ (we) ^ Pe imply that ^ a^ ^ 1 
for any J ^ J with Je > 0, and the inequality extends by continuity to the whole of J . We now 
calculate the derivative of a^ along Je for Je > and obtain, as 



d 



Pe 



that, for any J e J7with Je > 0, 

da-' 



P 



-^J,w 



n'^i^M) ^ti^e\Ay 



J,w 



dJe 

This implies in particular that 



= P< 1 



^r (^e) $r(^e|P7^) 



9a:^ 



9J. 



</3< 



and the comparison supj^^rp j^i Og ^ e^inf j__gro^ii Og follows. D 

Proof (Proposition 12 . 7p . According to Corollary 5.8 in [35] and to the Mean Value Theorem, we 



have 



Entp(exp(/A)) 
E(exp(A)) 



< 



4 ^ 

eeB(7?,) 



^ 



suP,/,e[o,i] 9^j exp(supj^g[o_i] fx] 
E(exp(/A)) 
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It is clear that 

dJ, 



— —Xpa^. 

•JJe 

On the other hand, Proposition 12.61 yields 

sup {aiY ^ sup flg s^ e'^ inf ai 

Jee[o,i] ./.e[04] Jee[o,i] 

and 

sup exp(/A) < e^^ inf exp(/A), 

J,e[04] -/.e[o,i] 

hence 

Entp(exp(A)) ^ A^/?^e^(i+^) ,p ^^ / .^^ ^,; ^, inf j.g[o,i] exp(A) 
E(exp(/A)) " 4 ^ Ue'lo.i]''^ "" E(exp(/A)) 

and the first bound follows. For the second one, remark that as we take infimums over Je we in 
fact obtain a quantity that is independent of Jg. Thus 

Entp(exp(/A)) A^^^e^(i+A) ,-^ /^ . ,, inf j^g[o,i] explA) 

E(exp(/A)) ^ 4 2^. Iv^^pj.'JIo,!]"^ "" E(exp(/A)) 

e£E{n) 



€ 



4top 




which ends the proof as 



\ee£;(7t) 

n 

Proof (Proposition 12. 8p . As 7?. = TZm = Ti-o. N,SN{S,n) is centered at the origin, we consider 
S^ = |o-:Z''^{±l} :cr^ = l,Vx^7^\a7^| 

the set of spin configurations on TZ with plus or mixed boundary conditions. The correspondence 
between the random-cluster representation (with q = 2) and Ising model gives 

T-' - 1op- •^ 

-n - ^,_i log ^,,± 

where Z^^ and Z.^^ are the partition functions 

0-GS+ \ e={x,y}e£(K) 

and Z.^"*" = ^ exp - ^ JeCfx^y , 

o-eS± \ e={a:,y}e£;(7t) / 

leading thus to 

(2.29) ai = n^^{(7^ay) - fi^ {<yxcry), Ve = {2:, y} E E{n) 

where /x.^^ is the Ising model on TZ with mixed boundary condition (plus on d~^TZ, minus on d~7V). 



We consider now an interface / for TZ as in Section 11.41 We recall that it is a minimal set of edges 
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such that connections from d^TZ to d TZ through E{TZ) \ I are impossible. We consider /+ the 
upper part of the interface /: 

1+ = {x : 3y e Z'^ : {x, y} e I and x ^ d^U in E{n) \ 1} 

and define symmetricahy the set I~ . We call then Sj the event that / is the spin interface between 
d~^TZ and d~TZ under the measure /i^ : 

,± cr(a:) = +l,Vxe /+ 
'^ ■ a{x) = -l,yx e I- 



Si = iae^i 



Conditionally on Si, the restriction of fi'^ to the upper (resp. lower) parts of TZ equals the Ising 
measure with uniform plus (resp. minus) boundary condition. Hence, for any {x, y} ^ I we have 

(2.30) fi^l (axCTylSi) ^ fj,T{^{aj;ay), 
and consequently 

(2.31) E «e ^ E 4^(5/)x2|/|. 

eeEilZ) I interface 

Thus it remains only to bound the average interface length, in the Ising sense, under /i^ . We 
remark that /i)^ {Si) can also be written as 

t^n {Si) - -j^ 

where ■^t?\/ stands for the partition function associated to the set of configurations with plus 

boundary condition on /+, I^ and on dTZ. Thanks to the assumption Jg ^ e and to the remarks 
that 

and Z^=^ ^ Z:l(+ cxp ('-/3|(9"7^|) , 

we have 

fi'^^iSi) s; exp(-/3£|/|+/3cdiV'^-i) 

as (5 < 1. We conclude with a Peierls estimate and bound the number of interfaces of cardinal 
n ^ 2cdN''-^/e by (c^)": 

^ al ^ ^^^^'^' + Y. n(crf)"e-'^-+''=''^'"' 

The second term goes to with N —^ 00 for f3 large enough. D 

Proof (Theorem II. 7p . The combination of Lemma 12.51 Propositions 12.71 and 12.81 implies that in 
the setting of Theorem 11.71 



< 



for P large enough, TZ ~ TZ^ — TZo,n,sn{S, n), S £ (0, 1) and N large enough. Integrating over A 
we obtain, as lim;^^o+ '''■r/^ = ^'^K' ^^"^ inequality 

C^/32g/3(l + A) 



Letting N ^ 00 gives 



'n ^ --'7^ -■ ^j 



C^/32 /3(1+A) 
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and the duality formula (|2.14|) yields the claim with c — CdJ"^"^/{(3'^ exp(2/3)), for large enough f3. 

a 

2.6. Concentration in a general setting. We give now the proof of Theorem 11.81 which is 
based on Herbst's argument, together with the controls of Lemma F2. 5 1 and Proposition 12.71 We 
will then give the proof of Corollarv ll.9l 

First we give an immediate consequence of the duality formula (I2.14p : 

Lemma 2.9. Assume that 

(2.32) limsup -7: ^ — c for some c ^ [0, 00]. 

Then, 

(2.33) lim sup -^"(^'(")-'^) ^ 1 g [0, 00] . 

Proof (Theorem II. 8p . Given (5 > and 5 G §„, we denote TZ^ the rectangular parallelepiped 
TZ'^ = TZo^n,sn{'>T',S) and introduce 

<^^ - liminf hminf 1 f^ ^^ G [0, 00] 

In view of Theorem 11.31 and Proposition 12. 21 we have 

T^(n)-AT«(n) == lim T^iv - AErijv 

Km A / ' A (ikL] dX' 



'0 

.J _ 



N^oo ,L dX' \ X' 



as Er.j^jv — lim;y_>o+ tj^n/X for any N finite. Lemma 12751 and Proposition 12 . 71 vield. for any e > 0: 

^X(^n)-XT^{n) ^ /3V(i+^),. . ^,. . ^ [^ 9t^^ dX' 
lim sup -^ ^ lim mt lim mt - / „^ —— 

;^^o+ ^ 4mp \-*o+ N^oo X Jq op X' 

fi2 I3{l+e) 

4mp "^ 
and an immediate application of Lemma 12.91 gives, after the limit £ — > 0, the lower bound: 

(2-34) hmsup ^ 

The lower bound is positive when K'^'^ < cx). In order to show that this is the case for Lebesgue 
almost all /3, we evaluate the integral of K^'^ on some interval [(3i, f32]- For any 6 > and 5 G Sn, 
Fatou's Lemma and Fubini Theorem imply that 

/"' Kl'f'dP ^ hm inf hm inf i /' /^ ^ ^ 

= liminfliminfi r^^^^^^^^dA'. 



A^0+ N^OO X Jq A' 

The convergence as TV — > 00 is uniformly dominated (by Jensen's inequality and Proposition 1 1.2[ 
^ T^ ^ Xcd/3) hence we finally obtain 

r^K^^^dp ^ liminf i/'*lz*), A' 
7/31 ^-0+ ^7o A' 

(2.35) = fl{n)-fl{n). 

in view of (|2.18|1 . In particular, K^'^ is finite for Lebesgue almost all /3 ^ 0. D 



26 MARC WOUTS 



We would like to make a remark on K^''^ . In view of Corollary II. 91 for Lebesgue almost every 
/3i, /32 with /3i < (32 one can replace fl (n)— fl (n) in (|2.35p with rl {n)—Trj^ (n). As a consequence, 
whenever T%{n) is derivable on some interval, K^'^ ^ dT'i{n)/d(3 for Lebesgue almost every /3 in 
that interval. 
Proof (Corollary II. 9p . We denote by 

T«_(n)=^hm T«_^(n) 

the left limit of T%{n). For any r G R, /3 i~* In{T) is non-decreasing hence flln) (defined at (|1.47p ) 
does not decrease with (i. According to Theorem II. 81 f%{n) coincides with rlln) for almost all /3, 
hence 

Tl-{n)^Tl{n)^Tl{n)yp^G 

hence the left continuity of rlln) at a particular (3 implies that f%{n) = t%{ti), in other words 
that lower deviations are (at least) of surface order. This is the first part of the claim. Now we 
consider 

and prove that T> is at most countable. The homogeneous extension of r? (n) to R'^ is convex as 
the pointwise limit of the /n_g, hence r? (n) is a continuous function of n e S"^^^. Consequently, 
for any dense sequence (wrOneN in 3"^^^, we have 



V C U {/3€R+:r|(n„)^r|_(n„)} 



nSN 

which is at most countable. D 

3. Low TEMPERATURES ASYMPTOTICS 

Here we study the low temperature asymptotics of surface tension and prove the results presented 
in Section 11.41 We begin with upper bounds on surface tension which hold in all generality, and 
then establish lower bounds with the help of Peierls arguments. 

3.1. Upper bounds on surface tension. Relevant upper bounds on surface tension are easily 
established: 

Lemma 3.1. Let F be a product measure P on [0, 1]"^, n E S"^^^ and A > 0. Then, 

(3.1) T^^in) s; p^^in) 

(3.2) and T^{n) ^ ||n||ixlog- / 

'^ texp(— ApJe) 

Proof We begin with the proof of (|3.ip and consider a rectangular parallelepiped TZ. With the 
notations of Section fl"^ for all interface / G T{TZ), the DLR equation yields 

<^'^- {Vn) ^ <^t iZi) ^ n *S(^'= = 0) = c^P {-PY.A 

eel \ eel ) 

and consequently ri-ji^ PfJ'Ti' which implies p.ip taking 7?. = 7?.^ = T^o, n,sn{S, n) and N ^ oo. 
Similarly, in view of the definition ()2.15p we have 

\eel 

\I\ 1 

^ L'i-i °^Eexp(-A/?Je) 

which yields p.2p if we choose for / the interface of smallest cardinal in 1(7?.) , which has a cardinal 
approximately |jn||ii''^^. D 



SURFACE TENSION IN THE DILUTE ISING MODEL. THE WULFF CONSTRUCTION. 27 



3.2. Quenched surface tension and msLxinial flows. We present here the proof of Proposition 
11.111 which is based on a control of the length of the interface, using a Peierls argument, and then 
the proof of Proposition II . 12l which uses a renormalization argument. 
Proof (Proposition 1 1 . 1 ij) . Given a rectangular parallelepiped TZ, we have 

iei{TZ) iex{n) ee/ 

We decompose the sum according to the length of the interface: for any c > ||n,||i, 

IeI{TZ):\I\<cL''-^ 

(3.3) + E gl^le-^^-^^= 

/eI(K):|/|^cL''-l 

The first term is not larger than 

{c^qT'^"'" exp {-PL'-'iii) 
and the expectation of the second one is 



^ E .l^'e^^^--^^ < i_.J.r.^...^ X ^^^^(^ 



.IeI(n):\I\'^cL'i- 






if p/3 = CrfgE(e '^'^=) < 1, which is the case for (3 large as P(Je = 0) = < (c^q) ^ . For any such 
/3, applying Markov's inequality we obtain, for any e > 0: 

^/GX(K):|/|^cL'i-i / ^^ 

Hence p.3p shows that, for J typical under P ~ up to large deviations of surface order - 

$^ {Vn) ^ icgr"^"'' exp (-/3L''- V^) + (P/^)'^"^)^^'" 
which proves that 

1 



T%{n) ^ min l3fi{n) - c\og{cdq),cAog 

for any /3 ^ such that pp < \. The lower bound is optimal for 

(3ii{n) 

C = T- 

log(cdg) + log — 

which is negligible with respect to (3 in the limit (3 -^ +oo, as log(l/p^) -^ +oo. The limit p.30p 
follows - the uniformity over n e S''*"^ is a consequence of the fact that /i is bounded. If Jmin > 0, 
then we even have, for some C < oo, that for f3 large enough (independent of n G S"^^^), 

T^in) ;^ Pp{n) - C. 

n 

Proof (Proposition I1.12p . The proof for (|1.3ip exploits a renormalization argument similar to 
the one used in [13]. As P(Je > 0) > Pc{d), for small enough e > it is still the case that 
P(Je ^ e) > Pc{d). We say that e £ E{Z'^) is open for J if Je ^ e, and consider the connected 
components for these definition of open edges. A block Bf = Ki + {1, . . . ,2i^}'' (i g Z'') of 
side-length 2K, is said good when, in Bf, 

i. there is a unique connected component of diameter larger or equal to K 
ii. and this connected component touches all the faces of the block. 
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The work of Pisztora [38] together with the knowledge that, in the case of independent bond 
percolation, the slab percolation threshold coincides with the threshold for percolation [55] imply 
that the sequence of random variables 

rji — Ijrpj^j, ^i^^y. gK ig good} 

stochastically dominate a site percolation process of parameter p close to one, provided that K is 
large enough. Provided that p (hence K) is large enough, a Peierls argument shows that there is 
a probability 1 — exp(cA^'^^^) (for some c > 0, for N large enough) that no X-block interface in 
TZ^ = TZq, n,sn{S, n) contains less than half of good blocks. 

Given a suitable K we now establish a lower bound on the quenched surface tension. Consider a 
realization J of the couplings with the property that no iiT-block interface in 7?.^ = Ti-o, n.sn{S, n) 
contains less than half of good blocks. For any such J, the event of disconnection requires the 
choice of a block surface of cardinality at least {N/KY~^, and that, in each good block, at least one 
edge with Jg ^ e be closed. Hence: for N large and J typical up to surface order large deviations. 



$^-(2?^«) ^ 



leading to t? ^ [(1 — e)l3e — log [c'^K'^q)] / K'^ ^ for large enough /?. The claim follows. 



D 



3.3. Surface tension under the averaged Gibbs measure. Under the assumption P(Je > 
0) = 1, we establish a lower bound on Tn[n) which is equivalent to the upper bound: 

Proposition 3.2. Assume that P(Je > 0) = 1. Then, uniformly over n £ S"^^^ , 
(3.4) T^{n) ^ (l-o^^eo(l))x ||n||i xlog 



Before we address the proof of Proposition 
consequence of Lemma 13.11 and Proposition 13.21 
Proof (Proposition 13. 2p . Remark that for any / C 1(7?.) , 



Eexp(-A/3Je) 
let us remark that Proposition 11.131 is a clear 



■/3Je 



ee/ 



ee/ 



If A < 1, the inequality (X]r=i ^0^ ^ SiLi ^\ foi" non-negative Xi yields 



$ 



7^ 



(^7^)'< E '^ni^i)\ ^ E ^'ni^it 



yiei{ii) 



hence 



E 



(*^"'(2?7^: 



< 



leiiii) 



Y: l[q'Ee~'P'' 
leiiK) eel 

/ei(K) 



As Ee~'^^'^= ^ as /3 ^- oo under the assumption P(Je > 0) = 1, the Peierls argument gives a 
relevant lower bound : there is c^ depending only on the dimension d, such that the number of 
interfaces of cardinal n in 1(7?.) is not larger than (c^)". Hence, 



E 



(^^"'(2?7^) 



€ 



€ 



n>minj-ez(7?,) l-f| 
1 



1 - Cdq^Ee-^PJ': 



[crfg^Ee-^'3j.]-i„,„<^,|/| 
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for /3 large enough, thus 



(3.5) 



T^^llnllix log 



1 



(1 + A)logCd- Alogg 



for all A ^ 1 and /3 large enough. If A ^ 1, Minkowski's inequality yields: 



E 



{'^i:"iT^n) 



1/A 



1/A 



^ ^ [£[($^-(2,) 

ieHTZ)eei 

and we conclude similarly that p.Sp holds again for all A ^ 1 and /3 large enough. The claim (|3.4p 
follows from the divergence lim^^+oo log(l/Ee~'^'''^'=) = +oo under the assumption P( Jg > 0) = 1, 
and the convergence is uniform in n e S"^"^ as p.Sp holds for any /3 large enough independent of 
n. D 



3.4. Limit shapes at low temperatures. The limit shape of Wulff crystals (|1.35p are 
immediately inferred from the uniform limits for surface tension. The Proposition below is a 
consequence of Proposition fTTTT] for the first point, Lemma [XT] and Proposition 13.21 for the second 



Proposition 3.3. Let P be a product measure on [0, 1]'' with P( Je > 0) = 1. 

i. The Wulff crystal W converges to the Wulff crystal W^ associated with the maximal flow ii. 
ii. For any A > 0, the Wulff crystal W^ converges to the hypercube Wll-lu = [±1/2]''. 

We remarked above that the maximal flow determines the limit shape of the crystal in the 
(quenched) dilute Ising model as the temperature goes to zero, while the crystals under the averaged 
Gibbs measure converges to the unit hypercube. 

Let us explain how the result of Durrett and Liggett [2T] for site first passage percolation can 
be used to show that VV^ is not in general an hypercube. We consider d — 2 and P such that 



(3.6) 



Je = ^ ) = P and P ( J, = 1) = 1 - p 



with pc < p < 1, where Pc is the critical threshold for oriented bond percolation. In the two 
dimensional case, one can also interpret /i(n) as the limit ratio over TV of the time needed for 
reaching the position Nn-^, if to every edge e we associate a passage time Je- Ii n belongs 
to the cone of oriented percolation (modulo the symmetries of the lattice 7?) one can reach the 
position Nn-^ following a directed path with all edges (except finitely many at the origin) satisfying 
Je = 1/2. Hence, in those directions. 



/i(n) 



1, 



" 1- 



However, the argument of [5T] (applied to bond in place of site first passage percolation) shows 
that when n is close enough to the axis, one has to use edges Je = 1 with a positive frequency, 
thus /x(n) > ||n||i/2 for those directions. The reciprocity formula 

T(n) = sup X ■ n 

for Wulff crystals, where W^ = {x : x ■ n ^ T{n),yn} is the un- normalized Wulff crystal for r, 
shows that W^ is not a square as /i(n) is not proportional to ||n||i. 
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4. Phase coexistence 

4.1. Profiles of bounded variation and surface energy. The coarse graining for the dilute 
Ising model (Theorem 5.10 in [33]) implies that at every position, the local magnetization Mk is 
close to ±771^ with large probability. In order to describe the geometrical structure of the phases, 
we estimate the probability that AI/c/to^ be close, in L^-distance, to a given Borel measurable 
function u : [0, 1]'' -^ {±1}- As a first step towards the description of phase coexistence, we define 
here the set of profiles we consider, define surface energy and the associated isoperimetric problem. 
In the following, £'^ stands for the Lebesgue measure on R'' and TC^"^ for the d—\ dimensional 
Hausdorff measure, which gives to any Borel set X C M"^ the weight 

W^-^ {X) = lim^ ^1^ inf \ YyAmm{Ei)Y-^ : supdiam(^,) <5,X ^\JeA 

where the infimum takes into account finite or countable coverings {Ei)i^i, and ad^i is the volume 
of the unit ball of R''"'^. The i^-distance between two Borel measurable functions u, v : [0, l]"* -^ M. 
is 



and the set L^ is 



\u-v\\li = / lu-vldC 

'[0,1]'* 



{u : [0, l]"* -^ R Borel measurable, ||u||ii < cx)} . 

In order that L^ be a Banach space for the L^-norm, we identify u : [0, l]"* ^ M with the class of 
functions v : |ju — w||ii — that coincide with m on a set of full measure. We also denote by V(u, 6) 
the neighborhood of radius ^ > in L^ around u G L^. 

For the study of phase coexistence, we have to consider virtually any u £ L^ taking values in 
{±1}. Before we can define the surface energy for such profiles, a description of the boundary of 
these profiles is necessary. It is done conveniently in the framework of bounded variation profiles 
(Chapter 3 in [4]). Given a Borel subset U C W^, the variation (or perimeter) of U is 

r{U) = sup Ij div/d/:^ / e cr (R^ [-1, i])| e [0, ex.] 

where C^(IR'^, [—1, 1]) is the set of C°° functions from W^ to [—1, 1] with compact support, and div 
the divergence operator: 

OXi OXn 

To U G W^ Borel measurable, we associate u = xu as in (|1.40p and define the set of bounded 
variation profiles BV as follows: 

BY = {u^Xu -U CI (0, 1)'' is a Borel set and V{U) < oo} . 

Bounded variations profiles u = xu G BV have a reduced boundary d*u and an outer normal 
n" : d*u -> S"*-! with, in particular, 'H'^-^{d*u) =V{U). 

This allows us to define the surface energy of bounded variation profiles. As the outer normal 
n" defined on d*u is Borel measurable, we can consider 



(4.1) T\u) 

and 



/ T'?«)d7i'^-i(x), VmgBV 



(4.2) J^^(w) = / T^{nl)d'H'^-^{x), Vu G BV,VA > 0. 

Jd*u 

where r'' (resp. r^) stands for the quenched surface tension of the dilute Ising model (resp. surface 
tension under the averaged Gibbs measure), see Theorem I 1 . 31 and (|1.24|) . Because the homogeneous 
extension of the surface tensions t* and t^ are convex (Proposition 12. 4|) , T"^ and T^ are lower 
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semi-continuous with respect to the L^-norm. See Chapter 14 in [TO] or Theorem 2.1 in [3J. For 
commodity, when u ~ xu G BV we also denote by !F'^{U) (resp. ^^{U)) the surface energy of u. 
When surface tension is positive, the level sets of T'^ and T^ are compact since, for all a ^ 0, 
the set 

(4.3) BVa - {u - xc/ e BV : r{U) ^ a} 

is itself compact for the i^-norm, cf. Theorem 3.23 in [4j. Consequently, T'^ and JF'^ are good rate 
functions. 

Let us conclude with a word on the solutions to the isoperimetric problem of finding the u G BV 
such that 

(4.4) / udC^i^— and J^9(w) is minimal ? 

The renormalized Wulff crystal VV (|1.35p is known to be the solution to the same problem without 
the constraint that U C (0, 1)''. Precisely, the solutions to [/ C M"^' Borel set with 

C^iU) == 1 and T'^iU) minimal 

are the translates of W', as the homogeneous extension of t'^ is convex (Proposition [2T4|) - see [40] . 
gH and [H]. 

For m < mp not too small, W' determines as well the optimal profiles in the cube (|4.4p . 
Consider a > with 

(4.5) "^ o 1 

The quantity a'' is precisely the least volume of U corresponding iou = xu & BV withj,p ^^^ udC^ ^ 

m/mp. If some translate of aW fits into the unit cube [0, 1]"^, that is if adiamoo(W^) ^ 1, then 
T{aW^) defined at p.4ip is not empty and therefore the infimum of !F''{u) for u G BV with 
/[o lid udC''^ ^ m/mp is exactly !F'^{ayV''). As a consequence, for all a satisfying adiamoo(W"?) ^ 1 

the optimal phase profiles correspond to the translates of aW that belong to [0, 1]'', which are 
the z + aW, for z G T{aW^). The same remains true if we replace ^f^ and W with T"^ and W^, 
for any A > 0. 

4.2. Covering theorems for BV profiles. Covering theorems play an essential role in the study 
of phase coexistence, as they allow to pass from the macroscopic scale (the phase profile u) to the 
microscopic scale (the dilute Ising model). We give first two definitions: 

Definition 4.1. Let u G BV, r : S"^^^ t-^ [0, oo] continuous, (5 > and TZ a rectangular parallelepiped 
as in (jl.lSp . included in [0, 1]''. We say that TZ is (5-adapted to u and r at x G d*u if the following 
holds: 
i. If n = n" is the outer normal to u at x, there are 5 G Sn and h G (0, S] such that, if 7?. C (0, 1)'' 
(we say that TZ is interior), then 

TZ ^x + hS+ [±Sh]n, 

and if 7?. n d[0, 1]'^ ^ (we say that TZ is on the border), then x G d[Q, 1]'^, n is also the outer 
normal to [0, l]"* at x and 

TZ^x + hS+ [Sh, 0]n. 
ii. We have 

H'^-^{d*uC^dTZ) =0, 

l--^rL'^~^{d*ur\TZ) s^5, 



and 

1 /■ 



^ / rin-)dTi'-^ 
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iii. If X : M'' ^ {il} is the characteristic function of 7?. defined by 



+ 1 ii{z-x)-n^O 
x\^) — "I 

then 



xi.) = \ :; ;;^^;--^"^^v..R^ 



'-.Jjx-um^^s. 



0. 



2Sh 

Definition 4.2. Let u £ BV, r : iS'^^^ i-^ [0, oo] continuous and (5 > 0. A finite sequence {Tli)i=i...n 
of disjoint rectangular paraUelepipeds included in [0, 1]'* is said to be a ^-covering for d*u and r if 
each TZi is (5- adapted to u and r and if 

(4.6) H^'-'Id^uXlJnA^s. 

The Vitali covering theorem (Theorem 13.3 in [10]) is especially well adapted to our purpose. 
Given a Borel set £■ C M'', we say that a collection of sets Z^ is a Vitali class for E if, for each 
X £ E and (5 > 0, there is U £ U with < diam U < 5 containing x. 

Theorem 4.3. [Vitali] Let E C M.'^ be T-C^^^ -measurable and consider U a Vitali class of closed 
sets for E. Then, there is a countable disjoint sequence (Ui)i^i in lA such that 

either ^(diamt/^)''"^ = oo or W^-^ ( ^ \ IJ ^M 

The Vitali Theorem allows us to state a short proof of the following: 

Theorem 4.4. For any u G BV, r : S'^^^ i-^ [0,oo] continuous and S,h> 0, there is a 6-covering 
{TZi)i=i...n for d*u and r. 

Before we give the proof of Theorem 14.41 we recall a property of the reduced boundary (see 
Theorem 3.59 in [4]): 

Lemma 4.5. Let u e BV. For all x e d*u, for all 5 G (0, 1), all S £ §„" one has 

lim -1-n''-' (d*unn^j,MS,K)) = i- 

Proof (Theorem l4.4p . We design a set E that has zero W*" ^-measure and such that the collection 
of closed rectangular parallelepipeds 

Us — {TZ (5-adapted to u and r at x G d*u \ E} 

is a Vitali class for d*u \ (E). This is enough to prove the claim: thanks to the Vitali 
covering Theorem, this implies the existence of a countable disjoint sequence {TZiji^i of ^-adapted 
rectangular parallelepipeds with either 

^(diam7^,)''"^ =oo or W^-^ {d*u\\Jni\ =0. 
iei \ iei ) 

The first case is in contradiction with the inequalities Xlh^^^H!^^^ {d*uC\'R,i) ^ 1 — (5 and 
■^d-i ^g*y^ ^ QQ^ hence the second is realized and the Theorem is proved. 

We define the set E by its complement in d*u: d*u \ E is the set of all x G d*u such that, for 
all S G Sri" , the following holds: 

i. If X G d[0, 1]'', then n" is the outer normal to [0, l]'' at x. 
a. The set {h>0: W^'^ {d*u n dn^,^h,Sh{S, n^)) > 0} is at most countable. 

iii. lim,,^o+ 71^ W''"' (d*unn^^h,5h{S,0) - 1- 

iv. lim,,_o+ 7r^/a*„n7t..,,,,(5,rxs)^("")^^''"' = ^(O- 
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This definition for E implies that Us is a Vitah class of closed sets for d*u \ E. We conclude the 
proof of Theorem 14.41 showing that E has zero 7i'*~^-measure, and more precisely that each of 
conditions (i)-(v) is true for (at least) H''~^-almost all x £ d*u: 
i. This condition holds for all x € d*u because of the inclusion U C (0, 1)'' if m = XU: ci. Theorem 

3.59 in [4. 
ii. Since the volume of d*u is zero, (ii) holds for all x. 
iii. Condition (iii) holds for all x S d*u in view of Lemma WM 
iv. It is a consequence of the strong form of the Besicovitch derivation theorem (Theorem 5.52 in 

[3]) together with Lemma BTSl that condition (iv) holds for 7i'^~-'^-almost all x e d*u. 
V. Condition (v) holds for all x G d*u, cf. Theorem 3.59 in [4j. 

D 

4.3. Lower bound for phase coexistence. Here we establish lower bounds for the probability 
of phase coexistence. In view of the applications, in particular to the control of the dynamics [43j 
or Chapter 4 in [JS], we establish it for a large class of profiles, that include Wulff crystals and 
shapes with C^ boundary. 

Proposition 14.81 below relates the probability of an event of disconnection along the boundary 
of a given profile to the surface tension t'\ for a given realization of the media. In Proposition [4?9l 
we show that conditionally on this event of disconnection, phase coexistence has large probability. 
Then we state in Proposition l4.10l a lower bound on the probability of phase coexistence for both 
quenched and averaged measures. 

Given some region [/ C (0, 1)'', iV e N* and J > 0, we consider E^ ' the set of edges at distance 
at most NVdS from NdU: 

E^^^ == (e e E"'{AN),d{e,NdU) sC NVdsl 

(see Figure [5]) and call 

^N,s _{ ^ o '^ Vx e Ajv \ iV[/, 2/ e Aat n iV[/ with \ 

the event that disconnection occurs around dU. In order to be able to control the probability of 
'Djj' , we introduce the following definition: 

Definition 4.6. We say that a profile u — xu is regular if 
i. U is open and at positive distance from the boundary d[0, 1]'' of the unit cube, 
ii. dU is d — 1 rectifiable and 
iii. for small enough r > 0, [0, 1]'' \ {dU + B{0, r)) has exactly two connected components. 

We recall that E C M.'^ is a. d — 1 rectifiable set if there exists a Lipschitzian function mapping 
some bounded subset of M.'^~^ onto E (Definition 3.2.14 in [23]). It is the case in particular of the 
boundary of non-empty Wulff crystals (Theorem 3.2.35 in [23]) and of bounded polyhedral sets. 
It follows from Proposition 3.62 in [4J that any u = xu regular belongs to BV and that dU — d*u 
up to a Ti*^^ ^-negligible set, so that the covering Theorem applies as well to dU . Assumption (ii) 
in Definition 14.61 has the following consequence: 



Lemma 4.7. Let u — xu G BV be a regular profile. Then, for any 5 > Q, for any 6-covering 
(711)1=1. ..n of u, one has 

Proof Clearly, the set 

n 

E = du\\Jn, 
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is a closed, d— 1 rectifiable set. Thus, the d— 1 Minkowski content of E equals the d— 1 dimensional 
Hausdorff measure of E (Theorem 3.2.39 in [53]). In other words: 

C'^{E + B{0,r)) 



lini 



2r 



n"^' (E) < 6 



and the claim follows. 



n 



Before we state Propositions 14.81 and 14.91 we give one more notation. The analysis of surface 
tension has been done for rectangular parallelepiped centered at lattice points. Changing the 
center of the parallelepipeds does not modify the behavior of surface tension, but this would have 
led to heavier notations. We prefer to proceed to a small adjustment here: given a macroscopic 
rectangular parallelepiped TZ C (0, 1)'' and N E N*, we let 

(4.7) 7^^ = 7V7^ + z^r(7^) 

where zn{TZ) G (—1/2,1/2]'* is chosen such that the center of VJ^ belongs to if. Still, for any 
finite collection (7?.i)i=i...„ of disjoint rectangular parallelepipeds in (0, l)"* and large enough N , 
the collection {TZf)i=i,,,n is disjoint and included in (0, A^)**. 



Proposition 4.8. Consider a regular u 
u, we have 



Xu ■ For any (5 > and any S-covering {TZi)i 



for 



(4.8) 



1 " 



cf3S 



for any N large enough, where c < oo depends on d and u. 



Proposition 4.9. Assume [3 > Pc o,nd (3 ^ J\f , and let u — xu regular. For any e > 0, for small 
enough 6 > there are K E N* and c > such that, for large enough N : 



(4.9) 



inf * 



J,w,+ 






M 



K 



mp 



£V{u,e) 



E. 



N.S 




^ e 



-cVN 



Proof (Proposition 



N,5 



To realize the event of disconnection T)^' , it is enough to realize all 



the Vj^N and to close all the edges that are at distance less than 1 + \/d from 

(n \ n 

du\\jTzAu\Jdutn, 
i=l / i=l 

where dutR- stands for the lateral boundary of TZ, that is the faces of &R. that are parallel 
to the orientation n of TZ. Thanks to Lemma 14.71 and Definition 14.11 there are at most 
ScdN"^^^ (l + 7i'*^^(9C/)) such edges for large enough N. An immediate application of the DLR 
equation yields (|4.8p . D 










/ 








1 / 













































Figure 5. The scales K and L 



N- 



Proof (Proposition 14. 9p . In order to obtain the claim for a mesoscopic scale K that does not 
depend on N, we proceed to a coarse grained analysis at two characteristic scales K and Ln = 
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[vN]. Given K e N*, we consider (A^, A^^g/^ ^ the {K,K)-coveTiiig of A^r as in Definition 5.f 
in [13] as well as the phase indicator 

given by Theorem 5.10 in [33], for the tolerance S. We call F — {0,1} '''^''^ the set of site 
configurations on the index of blocks Ikn,k- In order to apply the stochastic domination Theorem 
5.10 (iv) in [44j . we will define an increasing function f : F —^ {0, 1} with the appropriate properties. 
First, we need to describe the LAr-blocks: we call {Aj, A^)jgj„^ the (Ljv, LAr)-covering for Ikn,k 
as in Definition 5.1 in 0J. Then we let 

J = {j e Jnj< : Vz e K'^,E^ (A^) n E^^' = 0} 
and 

Given p G F a site configuration on Ikn.k and j £ J, we say that the L^r-block A' is good if 
there is a crossing cluster of open sites for p in A', of density at least 1 — 5. Then we define 
/ : F ^ {0, 1} letting 

i\P) — IjFor all j e J, A'- is good}' 

Clearly, / is an increasing function. We prove now that its expectation is close to 1 under high- 
parameter site percolation. Consider Bl, the site percolation process on / of density p £ (0, 1). 
According to Theorem 1.1 in [19 , for large enough p < I there is c > such that, for large enough 
N, for all j e J: 

B'p ({a; is good}) > 1 - exp {-2cL%-') 

and consequently (the cardinal of J is bounded by N'^), for p < 1 close enough to 1, for large 
enough N, 

Bl (/) ^ 1 - exp (~c^/N) . 

Consequently, the stochastic domination for (|(/>i|)ig/^ ^ (see Theorem 5.10 (iv) in [33]) yields the 
same lower bound on the expectation of /((|0i|)i£/): for large enough K (depending on 5), there 
is c > such that, for any N large enough: 



(4.10) Einf*;(-+J/((|0,|),,,,) 



iei ) 



^-cVN_ 



The event that / ((l^iDjgj) = 1 gives a control on the magnetization. For large enough N, the 
blocks (Ai)ig7 cover a fraction of Ajv that is close to 1 — C^ (dU + B(0, CdS)) — > 1. This and the 

properties of {(f>i)ieiA k (Theorem 5.10 (i) and (ii) in [33]) imply that, for small enough S > 0, for 
large enough N: 

f m\).ei) = 1 ^ ^ e Viu,e) OT^e V(l,e). 

We now consider a boundary condition vr e P^ ' . Because of the ^-disconnection, the spin of the 
clusters touching some A^ C NU with i £ I has a symmetric distribution under the conditional 
measure 



*A«,/3 



/ ((l^iDie/) = 1 and w = TT on E^'^j . 



Hence, one has 



inf ^^N.s ^«' 



./3(^eV(.,e) 



Ll) ~ n OH tjj 



^linf r^A«,4^((i'^'iW) 



cN.S 

uj — IT on Ejj 
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The claim follows as (|4l0| implies, as £^'' C E'^IAn) \ [j.^j E^^iA'^), that 



inf * 



J.+ 



, TreX) 



K,S ^N,P 



(fimhei) 



Lu — TT on £, 



N,5 



^ 1-e" 



c/2VlV 



sCe" 



C/2VJV 



D 



The final formulation of the lower bound for phase coexistence is the following: 



Proposition 4.10. Assume P > f3c and /3 ^ Af. For any ^ a < l/diamoo(W'') and e > there 
exists K (z fi* such that, 

(4.11) liminf — ^— -log^;('+ ( — - E V{xzo+aWi,£) ] > -■^''(XaWi) P-a.s. 

where zq = (1/2, . . . 1/2). Similarly, for any A > and ^ a < 1/ diamoo(W^), 



(4.12) 



1 



lim inf 



-logE 



y^N y^^^ixzo+c^w^^^] 



> -^\Xc.W^)- 



Proof Let U = zq + aW^. According to Theorem 3.2.35 in [53], dU is rectifiable, hence the 
profile u = xu is regular. Let e,S > 0. Thanks to Theorem 14.41 there exists a 5-covering {Tli)2^i 
adapted to the profile xu and r'. Proposition l4. 81 applies and gives, for (5 > small enough: 



J.+ i'Mk ^ ,,. 



--[/ 



(4.13) 



^ ^J^L*aV [-^^V{xv,e)\u = .on£-, 
xexp [-TV^-i [^/i^-1t^„ +c/3(5 j 



where c < oo depends on d and u. An important remark is that the two factors are independent 
under the product measure P. Proposition 14. 91 yields: 



(4.14) 



inf vI;;(-'+(':^ 



G V{xu, e)\u = TT on 8^^'^ < ^ U e'^^ . 



We prove first (|4.1ip and consider 7,^>0. If(5>0is small enough. Theorem 11.41 tells that the 
P-probability that t£jv > ''"'(wj) + 7 for some i e {1, . . . , n} decays like exp(— ciV*) where c > 0. 

Hence, with P-probability at least 1 — g^'^^^/'^ we have 



Nd- 



i=l 



for small enough 5 > and 7 > 0. Borel-Cantelli Lemma ensures that P-almost surely, 
liminf ^TT-YlogA*A^ ( -— ^ G V{xzo+aW-,£) ] > -^''iXaw) -C 



AT^oo N^' 



mp 



and (|4.1ip follows letting ^ ^ 0+. We conclude with the proof of (|4.12p . take A > and denote 
here U = zo + aW'^. Again, there exists a (5-covering {TZi)f^i adapted to the profile xu and r^. For 
N large enough, the 7?.f are disjoint and hence the t.£jv are independent under P. Consequently, 
for N large enough and A > 0, (liTT^)) and (|i?TIl) give 



E 



"i'^l^-^'--^ 



5s 



1 



2x3 



A 



J^Ecxp (-ATV^-i/if-V^ 



Kf 



X exp {~\N'^-^c(3S) 
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In view of Proposition 12. 2) this means 
1 



lim inf ■ 



-logE 



and the claim follows as S 



^i:(^^nxu,e) 



\ ^13 



> -^/i^V(n,)- Ac/3,5 



0. 



n 



4.4. Upper bound for phase coexistence. Here we address the opposite problem of providing 
an upper bound on the probability of phase coexistence along a given phase profile. Our analysis 
follows the same line as [S1E1[TT]. The cost of phase coexistence is easily related f Proposition 14. ll]l 
to another notion of surface tension (|4.15p . that uses a L^-characterization of phase coexistence. 
Then the L^-notion of surface tension is related to a percolative definition of surface tension with 
/ree boundary conditions, with the help of the minimal section argument f Proposition 14. 1^ . As 
in the uniform setting [llj . the surface tension with free boundary condition differs very slightly 
from the usual notion of surface tension f Proposition 14. 1 3p . 

The L^-dcfinition of surface tension is as follows. Given 6 > 0, a rectangular parallelepiped 
7^ C [0, l]'^ as in Definition [H] (i) and K,N eN* we define 



(4.15) 



~J,S.K 

'''ntz 



Wr^'^'Mll ''- 



NTZ 



M 



K 



mp 



^ 2SC'^ (7^) 



Li(K) 



where x is the characteristic function of TZ as in Definition 14.11 (iii) , and /i'^ the Gibbs measure 
on NTZ with boundary condition a. We have: 

Proposition 4.11. Let u G BV, 5 > Q and assume that {'R,i)i=i...n is a 5-covering for u. Then, 
for any e > small enough, any K, A^ G N* one has: 



(4.16) 



1 



irrlog/^A 



J,+ 



[M 



K 



J,,-, ---A« V ^^ e nu.e) ) < -Y.'^-'rii^ 



i=l 



Proof For e > small enough, the implication 
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K 
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G V(u,£) 
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K 
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i^SC^in,), ViG{l,...,n} 



ii(K,) 



holds. Thanks to (iii) in Definition 14.11 for such e we have 



M 



K 



mp 



G V(u,£) 



M 



K 



mp 



sC 2(5£'*(7^,), ViG {l,...,n}. 



ii(K.) 



J.+ 



Now, the Gibbs property for /i^' implies that 



/^A« 



M 



K 



mp 



&V{u,e] 



€ 




sC2J£''(7^,),ViG {l,...,n} 



s^ 2SC'^(n,) 



L^(1l, 



D 



^ exp 

thanks to ()4.15p . and the claim is proved 

Using the minimal section argument as in 5 one can compare the L^-surface tension to the 
surface tension under free boundary condition in 7?. = TZx.L.H{S,n), defined as 

(4.17) f^ = --^log$^^(I?^) 

where TZ = TZx^L,H/2{S,'n) is a rectangular parallelepiped twice finer than TZ. 
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Proposition 4.12. Assume f3 > (3c with (3 ^ TV. Then, there exists Cd.s G (0,oo) with 
lim^^o Cd,5 = such that, for any TZ as in Definition \4.1\ (i), for any 5 > 0, if K is large enough 
then: 



(4.18) limsup ^ logP (f^'^^ s; f^„ - cd,s) < 0. 



We do not detail here the proof of Proposition 14.121 as it is easily adapted from [5] . Then, the 
argument of [llj let us quantify the influence of the boundary condition on the value of surface 
tension: 

Proposition 4.13. Assume (3 > (3c o,nd (3 ^ M . Let TZ be a rectangular parallelepiped TZ as in 
Definition \4.1\ (i), with 6 G (0, 1). Then, 

(4.19) limsup —. logP (f^„ =^ r^„ - Cd5) < 

where Cd < oo depends on d only. 

We cannot afford to give here the proof of Proposition l4.13l as the generalization to the random 
case of the argument of |11| makes it far too long. However, no new ingredient needs to be 
introduced with respect to the original construction jll] , and the interested reader can consult the 
PhD thesis [3F for a complete development of the proofs of both Propositions 14.121 and 14.131 

The consequence of the three last Propositions, together with Varadhan's Lemma, is a lower 
bound on the probability of phase coexistence along a given profile under quenched and averaged 
measures; 

Proposition 4.14. For all (3 > (3c with (3 ^ M , for every u e BV and ^, A > 0, there exists e > 
such that, for K e N* large enough, 

(4.20) limsup-^ log A*:('; f:^ e V{u,e)] ^ -T^{u)+S, 



in ¥-probability (and P-almost surely if (3 ^ Afj) and 

1 A 



(4.21) linisup-^^YloglE 



N 



Nd 



' ■ ^'''^(t^'^ >'<«-)" ' -' 



^-T^{u)+C 



Proof We fix (5 S (0, 1) and a (5-covering {TZi)i=i,,,n for u as in Definition 14.21 We examine first 
the quenched convergence: according to Propositions 14. 12l and BT51 there is c > such that 

(4.22) P (f^'^^ > T^« - cd,s - cd<5) ^ 1 - exp(-ciV'^), Vi = 1 . . . n 

for K and N large enough. On the other hand, for any e > small enough Propositions HTTT] yields 

T log /^Aiv --— e V(m, e) < - 2^ K T^ii^ 



and hence, for K and N large enough, 

^logf.'^^ (^ e V{u,e)) ^-Y.f't' K^ - Cdj-CdS] 



with P-probability greater than 1 — nexp(— ciV''). This implies (I4.20p for (5 > small enough 
in view of the convergence t^n —^ T'^{ni) in P-probability (Theorem II. 3p or of the almost-sure 

convergence if /3 ^ Afi ( Corollary 1 1.9p . We examine now the averaged convergence: consider A > 
and again, a 5-covering {TZi)i^i,,,n for u. For K, N large enough and e > small enough we have 
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^ Ecxp -^A(/i,7V)^-ifi^f^ 



' NUi 
i=l 



< n eM-cN'^) + E exp - ^ Xih.Nf-^T^j. 



i=l 



e^p i Xj^ht^N''-^ icd,5 + CdS) 



X 

V 4=1 / 

in view of (|4.22p . Varadhan's Lemma f Proposition 12. 2p yields: for any e > small enough, any K 
large enough, 

r 1 1 w f I" "',+ f^K ^"'^ 

lim sup ^^ , , log ? ' ■ ■ ' ' 



<M^^^("'^) 



I 



<-E'^'"'[^^"0-Cd,.--Cd5] 



i=l 



and the conclusion follows for 5 > small enough. D 

4.5. Exponential tightness. The last step towards the proofs of Theorems II. 151 and II. 161 is the 
exponential tightness property. Note that the compact set BVa was defined at (14. 3p . 



Proposition 4.15. For any [3 > Pc with j3 ^ M , there exists C > and for every (5 > 0, Jar any 
K £W' large enough one has 

(4.23) hm^sup ^ log^liit. [^ i V(BV,, sA <k -Ca. 

The proof of Bodineau, loffe and Velenik given in [6] applies as well in the present case. 

4.6. Proofs of Theorems 11.151 to 11.201 Theorems 11.151 and 11.161 are consequences of the 
large deviations estimates (Propositions 14.101 and I4.14p together with the exponential tightness 
(Proposition I4.15P in view of the compactness of B Vq . The case of averaged Gibbs measures 
f Theorems II. 171 [1.191 and II. 20p presents complete similarity with the non-random case and for this 
reason we focus here only on the quenched case. Furthermore, the proof of Theorem 1 1.1 51 is similar 
to that of Theorem I1.16[ which is the reason for which we give the proof of p.42p only. 
Proof (First half of Theorem I1.16P . First we establish the lower bound 

(4.24) linnnf ^^ \og^i'f(^ (^ ^ 1 - 2aA ^ -J^'iXaW^), P-almost surely. 
The proof goes as follows: for any a' > a, for small enough e > one has 

^eV(x.o+.'w.,s)^^^^l-2a'^ 

hence, Proposition 14. 141 gives: for any a' > a, 

liminf j^ logA^A'^ (^ ^ 1 - 2aA ^ -^'(xa'wO, F-almost surely. 

The lower bound (|4.24p follows if we let a' -^ a. 

Now we establish the following upper bound: for any e > 0, there is (5 > such that 

^^^'^^P7Jd^^°Sl^AM — - i U ">^(X:r+aW<j,£) and ^^^^ =s: 1 - 2a'' 



(4.25) < -J'^iXa 



Y\>i , 
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in P-probability (P-almost surely if /3 <^ A//). To begin with, we choose a > so large that Ca in 
Proposition 14. 151 is larger than '2,T''{xaV\!i) + 2. Thanks to Markov's inequality, this implies that, 
for any 7 > 0, for large enough K^ 

(4.26) lim^sup ^ \ogii-'^+ (^^ i V(BV„ 7)) ^ -T'^iXo.w^) ~ 1, 

P-almost surely (see (14. 3p for the definition of BVa). Consider r/ > and let 



F = i u e BVa : / w s; 1 - 2a'^ + r/ and u ^ U "^ ( 



Xx+aWi, - 



For 7 > small enough, for large enough N the event 



implies that 



^ ^ U V(x.+aw., e) and ^^^ < 1 - 2a'' 



^^' -V{BY.n)or^eViF,j). 



The probability of the first event is under control (|4.26p for any 7 > (and large enough K) , hence 
we focus on the probability of the second one. Given ^ > 0, applying Proposition 14.141 we obtain 
e : u e BV i-^ e(u) G (0, S,) such that, for any u e BV and any K large enough: 

(4.27) hmsup -^^ logA.;('; (^ e V{u,e{u))] ^ -T^{u) + C 



N- 



Nd-i ^^^" V "1/3 



in P-probability (P-almost surely if /3 ^ A//). The set BVo is compact for the L^-norm, thus it 
can be covered by a finite union BVa C lj"=i V(ui,£(wi)) with Ui G BVq, i = 1. . .n. Since the 
right-hand side term is open, for 7 > small enough we still have 

n 

V(BVa,7)c [jV{u„e{u,)). 

i=l 

We consider (M'j)i=i...i the subsequence of the Ui such that V(ui,e(Mi)) intersects V(i^, 7). Thanks 
to the inclusion 

I 

V(F,7)cUv(uU«)) 
1=1 

and to (|4.27p . we have: for small enough 7, for large enough K: 

limsup — ^ log/i^ + I ^ e V(F, 7) ) < - inf T'^iu) + ^ 

in P-probability (P-almost surely if /? ^ A//). Yet, the limit as ^ ^ of the right-hand side is 
bounded from above by — inf„gi?/ T'^{u) where 



< u e BY a ■■ u ^ I - 2a'' + 2r] und u (^ I J V (xx^ 



for any 77 > 0. Yet, — inf^gi?' T'^{u) is strictly smaller, in the limit 77 — > 0, than -~T''{xa'wi) since 
the solutions to the isoperimetric problem (|4.4p are excluded. Together with (|4.26p . this implies 
(|i^ and the conclusion p"i^ follows from (ji?^ and (|i^ . D 
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4.7. Localization of the Wulff crystal under averaged measures. One consequence of the 
introduction of the random media is the localization of the WulfF crystal if the volume constraint 
acts on the media as well: the surface tension appears to be reduced on the contour of the crystal. 
Here we give the proof of Theorem 11.211 after a we state the following immediate consequence of 
the lower large deviations described in Theorem 11.61 



Lemma 4.16. Let TZ 



N 



Ti-o.N.SNiS^n) and j > 0, A 



(n) -7,'r^'+(n) +7]. Then, 



limsup— ^3Yl°gI^ 



N 



Nd 



\rl,eA^^}''''''Pi-^^'''^i-) 



<T^{n). 



Proof (Theorem ll.2ip . According to Theorems 11.171 and 11.201 it is enough to prove that 



limsup-^-^^3Ylog^ 



Af— >c» 
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J.+ 



l-flN 



e A" and 



M 



K 



mp 



Xz+aW^ 



^e 



(4.28) 



< -T^{aW^ 



In the case that the parallelepiped TZ. does not intersect the crystal z + adW'^, for 5 > small 
enough any ^-covering {TZi)i=i,,,n for z + aW^ and r"? does not intersect TZ. For e > small 
enough and K large enough, Propositions 14.111 14.121 and 14.131 the definition of the (S-covering and 
the independence of t^ from the Tj^-j^. under the product measure P imply that the right-hand 
side of (|4.28p is bounded from above by 

-^\aW^) + o (1) + hmsup -i^ logP (r^„ e ^^) 



N 



Nd- 



which is strictly smaller than —!F^{aW''^) for small enough 5, as the last term is strictly negative. 
Now we consider the case when the parallelepiped TZ is tangent to the crystal. For h > small 
enough, for e > small enough, the strict inequality 



lim sup ^^^_^ log E 
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Nd- 



sup fl 



,J,a 
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mp 



Xz+aW>- 



^ e 



L^{n) 



(4.29) 



< - / T^{n,)dU 



holds according to Propositions 14. 12l and l4T3l and Lemma F4. 161 Let {TZi)i=i,,,n be a 77-covering 
for z + aW^ and r^. Propositions 14. 11"] 14. 121 and 14. 131 and the properties of the 7y-covering imply 
that for £ > small enough (depending on 77) and large enough K, the cost of phase coexistence 
outside of TZ is bounded above by 



limsup^^Tj^Ylogl^ 
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Nd- 
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sup /i 



vi:K,nK=0'^est^ 



,J,S 



M 



K 



mp 
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s; e 



mn) 



< 



T^in,)dn 



d{z+aW^)\K 






Thus, choosing h > small enough then e > small enough and K large enough, the strict 
inequality holds in (|4.28p and the claim follows. D 



42 MARC WOUTS 



References 



[1] M. Aizenman, J. T. Chayes, L. Chayes, and C. M. Newman. The phase boundary in dilute and random Ising 

and Potts ferromagnets. J. Phys. A, 20(5):L313-L318, 1987. 
[2] M. Aizenman, J. T. Chayes, L. Chayes, and C. M. Newman. Discontinuity of the magnetization in one- 
dimensional l/\x — y\'^ Ising and Potts models. J. Statist. Phys., 50(l-2):l-40, 1988. 
[3] L. Ambrosio and A. Braides. Functionals defined on partitions in sets of finite perimeter. II. Semicontinuity, 

relaxation and homogenization. J. Math. Pures Appl. (9), 69(3) :307— 333, 1990. 
[4] L. Ambrosio, N. Fusco, and D. Pallara. Functions of bounded variation and free discontinuity problems. Oxford 

Mathematical Monographs. The Clarendon Press Oxford University Press, New York, 2000. 
[5] T. Bodineau. The Wulff construction in three and more dimensions. Comm. Math. Phys., 207(l):197-229, 1999. 
[6] T. Bodineau, D. loffe, and Y. Velenik. Rigorous probabilistic analysis of equilibrium crystal shapes. J. Math. 

Phys., 41 (3): 1033-1098, 2000. Probabilistic techniques in equilibrium and nonequilibrium statistical physics. 
[7] B. BoUobas. Graph theory, volume 63 of Graduate Texts in Mathematics. Springer- Verlag, New York, 1979. An 

introductory course. 
[8] P. Carmona, F. Guerra, Y. Hu, and O. Menjane. Strong disorder for a certain class of directed polymers in a 

random environment. J. Theoret. Probab., 19(1):134-151, 2006. 
[9] R. Cerf. Large deviations for three dimensional supercritical percolation. Asterisque, (267):vi-|-177, 2000. 
[10] R. Cerf. The Wulff crystal in Ising and percolation models, volume 1878 of Lecture Notes in Mathematics. 

Springer- Verlag, Berlin, 2006. 
[11] R. Cerf and A. Pisztora. On the Wulff crystal in the Ising model. Ann. Probab., 28{3):947-1017, 2000. 
[12] D. Chafa'i. Inegalites de Poincare et de Gross pour les mcsurcs de Bernoulli, de Poisson, et de Gauss. 

Unpublished, http: //hal. archives-ouvertes. fr/ccsd-000124.28 2005. 
[13] J. T. Chayes and L. Chayes. Bulk transport properties and exponent inequalities for random resistor and flow 

networks. Comm. Math. Phys., 105(1):133-152, 1986. 
[14] J. T. Chayes, L. Chayes, and J. Frohlich. The low-temperature behavior of disordered magnets. Comm. Math. 

Phys., 100(3):399-437, 1985. 
[15] L. Chayes, J. Machta, and O. Redner. Graphical representations for Ising systems in external flelds. J. Statist. 

Phys., 93(l-2):17-32, 1998. 
[16] F. Comets, T. Shiga, and N. Yoshida. Probabilistic analysis of directed polymers in a random environment: a 

review. In Stochastic analysis on large scale interacting systems, volume 39 of Adv. Stud. Pure Math., pages 

115-142. Math. Soc. Japan, Tokyo, 2004. 
[17] F. Comets and V. Vargas. Majorizing multiplicative cascades for directed polymers in random media. ALEA 

Lat. Am. J. Probab. Math. Stat., 2:267-277 (electronic), 2006. 
[18] A. Dembo and O. Zeitouni. Large deviations techniques and applications, volume 38 of Applications of 

Mathematics (New York). Springer- Verlag, New York, second edition, 1998. 
[19] J.-D. Deuschel and A. Pisztora. Surface order large deviations for high-density percolation. Probab. Theory 

Related Fields, 104(4):467-482, 1996. 
[20] R. Dobrushin, R,. Kotecky, and S. Shlosman. Wulff construction: A global shape from local interaction, volume 

104 of Translations of Mathematical Monographs. American Mathematical Society, Providence, RI, 1992. 
[21] R. Durrett and T. M. Liggett. The shape of the limit set in Richardson's growth model. Ann. Probab., 9(2):186- 

193, 1981. 
[22] R. G. Edwards and A. D. Sokal. Generalization of the Fortuin-Kasteleyn-Swendsen-Wang representation and 

Monte Carlo algorithm. Phys. Rev. D (3), 38(6):2009-2012, 1988. 
[23] H. Federer. Geometric measure theory. Die Grundlehren der mathematischen Wissenschaften, Band 153. 

Springer- Verlag, New York, 1969. 
[24] I. Fonseca. The Wulff theorem revisited. Proc. Roy. Soc. London Ser. A, 432(1884):125-145, 1991. 
[25] I. Fonseca and S. Miiller. A uniqueness proof for the Wulff theorem. Proc. Roy. Soc. Edinburgh Sect. A, 

119(1-2):125-136, 1991. 
[26] G. Grimmett and J. M. Marstrand. The supercritical phase of percolation is well behaved. Proc. Roy. Soc. 

London Ser. A, 430(1879):439-457, 1990. 
[27] D. A. Huse and C. L. H. Henley. Pinning and roughening of domain walls in Ising systems due to random 

impurities. Phys. Rev. Let., 54(25):2708-2711, 1985. 
[28] D. loffe. Large deviations for the 2D Ising model: a lower bound without cluster expansions. J. Statist. Phys., 

74(l-2):411-432, 1994. 
[29] D. loffe. Exact large deviation bounds up to Tc for the Ising model in two dimensions. Probab. Theory Related 

Fields, 102(3):313-330, 1995. 
[30] D. loffe and R. H. Schonmann. Dobrushin-Kotecky-Shlosman theorem up to the critical temperature. Comm. 

Math. Phys., 199(1):117-167, 1998. 
[31] K. Johansson. Shape fluctuations and random matrices. Comm. Math. Phys., 209(2):437-476, 2000. 
[32] H. Kesten. Aspects of flrst passage percolation. In Ecole d'ete de probabilites de Saint-Flour, XIV — 1984, 

volume 1180 of Lecture Notes in Math., pages 125-264. Springer, Berlin, 1986. 



[33 

[34 

[35 

[36 

[37 

[38; 

[39 

[40" 
[41 

[42; 
[43; 

[44 

[45; 



[46' 

[47; 



SURFACE TENSION IN THE DILUTE ISING MODEL. THE WULFF CONSTRUCTION. 43 



H. Kesten. Surfaces with minimal random weights and maximal flows: a higher-dimensional version of first- 
passage percolation. Illinois J. Math., 31(1):99— 166, 1987. 

J. F. C. Kingman. The ergodic theory of subadditive stochastic processes. J. Roy. Statist. Soc. Ser. B, 30:499- 
510, 1968. 

M. Ledoux. The concentration of measure phenomenon, volume 89 of Mathematical Surveys and Monographs. 
American Mathematical Society, Providence, RI, 2001. 

A. Messager, S. Miracle-Sole, and J. Ruiz. Convexity properties of the surface tension and equilibrium crystals. 
J. Statist. Phys., 67(3-4) :449-470, 1992. 

C.-E. Pfister. Large deviations and phase separation in the two-dimensional Ising model. Helv. Phys. Acta, 
64(7):953-1054, 1991. 

A. Pisztora. Surface order large deviations for Ising, Potts and percolation models. Probab. Theory Related 
Fields, 104(4):427-466, 1996. 

R. Rossignol and M. Theret. Lower large deviations for maximal flows through a box in first passage percolation. 
arXiv:0801.0967, 2008. 

J. E. Taylor. Crystalline variational problems. Bull. Amer. Math. Soc, 84(4):568-588, 1978. 
M. Theret. On the small maximal flows in first passage percolation. To appear in Annales de la faculte des 
sciences de Toulouse, 2007. 

M. Theret. Upper large deviations for the maximal fiow in first passage percolation. Stochastic Process. AppL, 
117{9):1208-1233, 2007. 

M. Wouts. Glauber dynamics in the dilute Ising model below Tc. In preparation. 

M. Wouts. A coarse graining for the Fortuin-Kasteleyn measure in random media. Stochastic Process. AppL, 
doi:10.1016/j.spa.2007.11.009, 2007. 

M. Wouts. The dilute Ising model : phase coexistence at equilibrium & dynamics in the region of phase 
transition. PhD thesis, Universite Paris 7 - Paris Diderot, |http : //tel . arch ives- ouvertes . fr/tel- 002728 99 [ 
2007. 

Y. Zhang. Critical behavior for maximal flows on the cubic lattice. J. Statist. Phys., 98(3-4):799— 811, 2000. 
Y. Zhang. Limit theorems for maximum flows on a lattice. Arxiv:0710.4589, 2007. 



Modal'X, Universite Paris Quest - Nanterre La Defense. Bat. G., 200 avenue de la Republique, 92001 
Nanterre Cedex. 

E-mail address: marc.woutsOu-parislO.fr 



